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ABSTRACT 


This  paper  presents  a  package  of  nine  programs  for  the 
TI-S9  calculator.  This  package  was  developed  as  a  solution 
to  two  problems.  One  problem  involved  expanding  and  modify¬ 
ing  an  existing  set  of  programs;  and  a  second  problem  involved 
developing  five  distribution  approximating  programs.  The 
solution  to  these  problems  represents  a  package  with  consider¬ 
able  capability  in  computing  confidence  intervals,  performing 
hypothesis  tests  and  approximating  distribution  values.  The 
distribution  approximations  include  inverse  CDF  values  for 
the  Normal,  Chi-square,  Student's  t  and  F  distributions, 
which  allow  the  computation  of  confidence  intervals  without 
using  tables. 

The  TI-59  proved  to  be  a  useful  tool  in  solving  these 
problems  and  demonstrated  the  capability  of  hand-held  programma¬ 
ble  calculators.  The  comprehensive  set  of  user  guides  included 
in  this  programming  package  provides  even  the  inexperienced 
user  with  a  step-by-step  introduction  to  this  capability. 
Additionally,  the  methods  used  in  preparing  this  programming 
package  are  directly  applicable  to  other  calculators  or 
computers . 
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I .  INTRODUCTION 


The  purpose  of  this  paper  is  to  trace  the  development  of 
a  package  of  nine  programs  for  use  in  the  TI-59  calculator. 

Two  of  these  programs  compute  confidence  intervals  for  either 
the  one -population  or  two -population  situations.  Another  two 
programs  perform  hypothesis  testing,  again,  for  either  one- 
or  two-population  situations.  The  remaining  five  programs 
generate  approximate  distribution  values  for  the  Normal, 
Binomial/Multinomial,  Chi-square,  Student's  t  and  F 
distributions . 

The  basis  for  this  programming  effort  was  a  set  of  four 
TI-59  programs  written  by  Professor  P.  W.  Zehna  for  his 
personal  use  and  later  used  in  the  classroom.  Professor 
Zehna 's  programs  also  computed  confidence  intervals  and 
performed  hypothesis  testing.  However,  his  programs  were 
not  completely  user-friendly,  especially  in  terms  of  user 
guides;  and  they  were  dependent  on  obtaining  some  percentile 
values  from  standard  tables.  This  paper  then  presents  a 
significant  expansion  in  the  scope  of  these  early  programs. 

The  main  thrust  of  this  expansion  includes  a  simplified  and 
standardised  set  of  programs  and  user  guides,  while  eliminating 
the  dependence  on  distribution  tables. 

This  package  of  nine  programs  was  designed  for  two  types 
of  users:  the  student,  who  might  be  asked  to  quickly  solve 
several  very  different  problems  in  succession;  and  the  working 
analyst,  whose  main  concern  is  one  specific  problem  which 


rx 
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requires  great  accuracy.  Of  particular  importance  to  both 
types  of  users  is  a  detailed  set  of  user  guides  (Appendix  A) 
which  includes  sample  problems.  Many  of  these  sample  problems 
have  been  solved  in  two  ways;  the  first  way  involves  using 
program-generated  percentiles;  the  other  way  requires  input 
of  tabled  percentiles.  While  it  is  not  absolutely  necessary 
to  look  up  tabled  values  when  using  the  confidence  interval 
programs,  that  option  has  been  included,  and  should  be  used 
when  increased  accuracy  is  desired.  Except  for  some  cases 
involving  small  degrees  of  freedom,  however,  the  accuracy  of 
these  first  two  programs  is  quite  good  as  can  be  seen  in  the 
sample  problems  of  Appendix  A.  When  tables  are  not  available, 
the  percentile  values  generated  by  the  five  distributional 
programs  is  outstanding  (Appendix  B)  and  can  be  used  in 
lieu  of  tabled  values.  One  very  convenient  feature  of  the 
distribution  programs  is  the  ability  to  provide  values  normally 
available  only  by  interpolation  in  the  standard  tables.  It 
should  be  noted  that  the  methods  used  to  generate  approxima¬ 
tions  in  the  confidence  interval  program  is  slightly  less 
accurate  than  those  used  in  the  distribution  programs.  This 
difference  is  due  to  the  limited  number  of  program  steps 
available  in  the  first  two  programs.  Except  for  the  two 
hypothesis  testing  programs,  there  is  no  requirement  to  use 
the  applied  statistics  module  in  the  TI-59.  However,  this 
module  is  required  when  performing  hypothesis  testing  because 
of  the  large  program  size  and  the  need  to  provide  a  signifi¬ 
cance  level  with  each  test.  Thus,  it  was  possible  to  provide 
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detailed,  yet  simple,  user  guides  to  implement  theoretically 
correct  and  accurate  programs.  By  providing  the  option  of 
using  either  approximations  or  tabled  values,  both  the  student 
and  the  analyst  can  accurately  solve  a  variety  of  problems. 


II.  NATURE  OF  THE  PROBLEM 


The  programming  effort  presented  in  this  paper  was 
generated  by  two  major  problems.  The  foremost  problem 
involved  expanding  and  modifying  an  existing  set  of  TI-59 
programs  into  a  more  user-friendly  package.  A  second, 
closely  related  problem  was  adding  the  capability  for  generating 
accurate  distribution  approximations.  The  following  discussion 
of  the  solutions  to  these  two  problems  is  a  general  overview  of 
the  particular  solutions  used.  The  specific  methods  and  theory 
of  solution  are  left  for  later  sections  in  this  paper. 

A.  USER-FRIENDLY  PROGRAMMING 

User-friendly  programming  implies  programming  with  the 
user's  knowledge,  ability,  and  hardware  familiarity  in  mind. 

In  this  light,  a  user-friendly  program  is  a  program  which 
has  significant  capability,  yet  can  be  used  by  those  with 
only  a  modest  knowledge  of  either  the  calculator  or  the  theory 
involved.  Applying  this  definition  as  a  framework  for  providing 
a  user-friendly  programming  package  resulted  in  four  areas  of 
effort.  These  areas  are:  standardized  data  entry,  standardized 
solution  procedures,  maximized  use  of  calculator  capabilities, 
and  improved  user  guides. 

1 .  Standardized  Data  Entry 

The  nine  programs  in  this  package  all  require  some 
form  of  data  entry.  The  data  entry  schemes  for  the  confidence 
interval  and  hypothesis  testing  programs  need  careful 
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standardization  because  of  the  similarity  of  the  data  and  the 
size  of  the  data  sets  involved.  There  are  three  possible 
types  of  data  in  these  first  four  programs.  Data  may  be 
from  a  one -population  sample,  a  two -population  paired  sample, 
or  a  two-population  independent  sample.  Any  of  these  data 
can  be  easily  entered  in  the  form  of  summary  statistics, 
if  available;  however,  raw  data  requires  some  standardization 
between  programs.  The  data  entry  schemes  of  the  first  four 
programs  use  similar  data  entry  subroutines  which  take 
advantage  of  the  I-59's  data  entry  sequence.  In  each  of 
these  programs,  the  data  entry  subroutine  is  initiated  by 
pressing  m.  The  sequence  of  data  entry  then  differs 
slightly  depending  on  the  type  of  data  being  entered.  All 
of  these  entry  subroutines  use  a  format  which  requires  a 
R/S  to  be  pressed  after  each  data  point  entry.  This  method 
saves  one  keystroke  for  each  data  point  compared  to  the  TI-59's 
two-keystroke,  2nd  Z  +  ,  entry  method.  Also,  the  [r/ sj 
key  is  very  close  to  the  numerical  keyboard,  compared  to  the 
2nd}  key,  thereby  eliminating  a  source  of  data  entry  errors 
that  frequently  occur  when  using  the  find}  113  sequence . 

The  data  entry  schemes  used  for  the  five  distribution 
programs  presented  no  standardization  problem.  All  of  these 
programs,  except  the  Multinomial,  require  only  a  few  parameter 
entries.  In  the  Multinomial  program,  which  can  accept  as  many 
as  35  pairs  of  parameters,  the  data  entry  problem  was  more 
difficult.  Each  of  the  multinomial  data  points  is  stored 
separately  until  computation  begins.  The  data  entry 
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sequence  automatically  repartitions  the  calculator  to  make 
room  for  this  potentially  large  amount  of  data. 

Data  stored  directly  in  registers  has  been  somewhat 
standardized  between  programs  by  using  the  same  registers, 
where  possible,  for  similar  data.  More  information  concerning 
the  contents  of  data  storage  registers  can  be  found  in 
Appendix  A. 

2 .  Standardized  Solution  Procedures 

A  considerable  effort  was  made  to  standardize  the 
steps  required  in  each  problem  solution  sequence.  In  the  two 
confidence  interval  programs  (see  Appendix  A)  problems  are 
solved  in  three  steps.  In  the  first  step,  the  calculator  is 
repartitioned  and  the  data  are  entered.  The  second  step 
requires  entry  of  either  a  percentage ,  in  which  case  the 
program  generates  an  approximate  percentile,  or  a  previously 
obtained  percentile  value.  The  last  step  involves  the  selec¬ 
tion  and  initiation  of  the  proper  solution  subroutine.  A 
similar  three-step  sequence  is  used  in  the  two  hypothesis 
testing  programs  where  the  first  step  includes  repartitioning 
and  entry  of  test  parameters.  The  second  step  involves  data 
entry,  and  the  third  step,  subroutine  initiation. 

The  five  distribution  programs  contain  only  a  limited 
amount  of  standardization  due  to  the  different  nature  of  the 
programs.  The  Normal,  Chi-square,  Student's  t,  and  F  programs 
use  the  EJ  DO-  □□  labels  for  the  same  basic  functions. 

The  ra  label  is  used  in  the  Normal,  Chi-square,  and  Student's  t 
(not  F)  to  generate  approximate  density  values.  The  [  B  ]  label 
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is  used  to  generate  CDF  approximations  in  all  four  of  these 
programs,  where  similarly,  the  [  C  f  label  is  used  to  generate 
inverse  CDF  approximations.  The  Binomial/Multinomial  program 
shares  none  of  these  standardized  label  uses. 

3 .  Maximized  Use  of  Calculator  Capabilities 

The  ability  of  the  TI-59  to  repartition  was  the  basis 

for  the  expansion  in  capability  over  that  achieved  by 

Professor  Zehna's  original  programs.  Repartitioning  allows 

the  addition  of  the  program  steps  necessary  for  the  inverse 

CDF  approximating  subroutines  in  the  confidence  interval 

programs.  The  added  programming  space  was  also  used  to  com- 

2 

pute  one  confidence  interval  estimate,  for  (J  with  JJ.  known, 
which  was  not  available  in  the  original  programs.  Additionally, 
repartitioning  makes  possible  an  F  distribution  program  and  a 
Multinomial  distribution  program. 

However,  repartitioning  is  not  without  its  price.  The 
larger  programs  now  require  three  edges  of  the  magnetic 
program  cards,  which  presents  a  slight  inconvenience  in 
added  loading  and  storage  requirements.  Other  problems 
associated  with  repartitioning,  such  as  inadvertent  loss  of 
program  steps  and  unwanted  or  improper  partitioning,  have 
hopefully  been  eliminated  from  this  programming  package  by 
extensive  validation  with  sample  problems. 

4 .  Improved  User  Guides 

Preparation  of  improved  user  guides  was  a  key  element 
in  making  a  user-friendly  programming  package.  The  complexity 
of  the  programs  involved  and  the  intended  use  by  students 
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necessitated  a  departure  from  the  standard  TI-59  program 
record  sheet.  As  can  be  seen  in  Appendix  A,  the  improved 
user  guides  are  organized  with  the  user  in  mind.  The  general 
outline  for  each  user  guide  follows  this  pattern: 

a.  Introduction 

b.  General  Procedures 

c.  Specific  Procedures 

d.  Additional  Capabilities 

e.  Labels  Used 

f.  Storage  Register  Contents 

g.  Sample  Problems 

The  user  guides  for  the  five  distrubiton  programs  have  been 
combined  to  take  advantage  of  the  relative  simplicity  in 
using  these  programs.  The  answers  to  each  sample  problems 
are  in  a  10-digit  format  to  provide  a  positive  check  on 
calculator  output  when  working  each  problem.  The  confidence 
interval  sample  problems  are  solved  in  two  ways  to  demonstrate 
the  differences  between  using  approximations  and  tabled 
values  for  the  required  percentiles. 

B.  ACCURATE  DISTRIBUTION  APPROXIMATIONS 

The  second  major  problem  addressed  in  this  section 
involves  generating  distribution  approximations  for  both  the 
confidence  interval  programs  and  the  distribution  programs. 

The  capability  to  compute  accurate  distribution  approximations 
provides  a  new  dimension  to  this  programming  package  by 
eliminating  the  need  for  standard  distribution  tables  when 
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solving  confidence  interval  problems.  A  drawback  of  this 
new  capability  is  the  time  required  to  obtain  some  values 
from  the  approximating  programs.  The  time  required  by  the 
confidence  interval  approximations  is  not  nearly  as  long  as 
that  for  the  distribution  programs;  but  then,  the  quality  of 
the  approximations  is  not  as  good  either  (see  Section  III). 

Appendix  B  contains  comparisons  of  tabled  values  with 
both  the  confidence  interval  program  approximations  (Type  I) 
and  the  distribution  program  approximations  (Type  II).  The 
actual  probability  values  used  in  these  comparisons  were 
obtained  using  the  distribution  programs  and  can  be  regarded 
as  being  very  close  to  the  actual  probability  achieved.  The 
missing  values  in  the  inverse  F  comparison  are  due  to  the 
inability  of  the  Type  I  approximation  to  generate  inverse  F 
values  when  either  of  the  degrees  of  freedom  parameters  is  one. 
Also,  while  only  selected  approximations  are  listed  in  the 
comparison,  the  inverse  CDF  approximations  are  of  nearly 
equal  quality  over  the  entire  range  of  the  appropriate 
function.  The  Type  I  approximations  displayed  in  Appendix  B 
are  generally  not  as  accurate  as  Type  II  approximations  in 
terms  of  actual  probability  achieved.  Only  the  inverse  CDF 
approximations  for  the  Chi-square,  t,  and  F  distributions  are 
presented  in  Appendix  B  since  all  the  other  approximations 
which  are  available  in  the  distribution  programs  duplicate 
table  entries.  These  other  approximations  include  probability 
values,  CDF  values,  and  various  other  distribution  values 
(see  Appendix  A) .  The  superior  quality  of  these  approximations 
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can  be  attributed  to  using  the  same  approximating  methods 
used  in  the  TI-59  Applied  Statistics  Module  (see  Section  III). 
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III.  THEORY 


This  section  on  theory  is  presented  as  a  background  for 
the  solution  methods  used  in  the  accompanying  programs.  As 
such,  it  is  not  intended  to  be  a  primer  in  statistics. 

Instead,  this  section  should  be  considered  as  an  intermediate 
level  derivation  of  the  specific  statistical  methods  used  in 
programming.  For  a  more  basic  explanation  of  this  material, 
the  references  cited  within  each  subject  area,  or  equivalent 
texts,  should  be  consulted. 

There  are  five  subject  areas  discussed  in  this  section. 
First,  the  theory  used  for  estimating  confidence  intervals 
in  the  first  two  programs  is  discussed.  Next,  the  hypothesis 
testing  theory  necessary  for  programs  three  and  four  is  dis¬ 
cussed.  And  lastly,  the  methods  used  in  all  nine  programs 
to  obtain  approximations  to  distribution  values  are  discussed. 

A.  THEORY  FOR  ONE  -  POPULATION  CONFIDENCE  INTERVAL  ESTIMATION 

The  derivation  of  theoretical  interval  estimates  will  be 
done  in  the  same  order  as  these  estimates  appear  in  the  User 
Guide  for  Program  1  (Appendix  A).  Most  of  these  derivations 
use  the  pivotal-qualtity  method  to  obtain  confidence  intervals 
(l,u)  [Ref.  1:  pp.  379-389],  Other  methods  used  here  will 
be  discussed  in  slightly  more  detail,  but  will  still  be  brief 
compared  to  the  referenced  texts.  When  forming  a  C.I.  with 
a  nonsymmetric  distribution  the  interval  will  represent  an 
equal  tails  solution,  where  equal  tails  implies 
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P [X  <  1]  =  P[X  >  u]  =  (1  -*Y)/2.  This  method  does  not  provide 
the  shortest  C.I.  for  a  given  7  ;  however,  this  method  is 
commonly  used  for  its  ease  of  computation  [Ref.  1:  p.  382]. 
Regardless  of  the  method  used,  the  resulting  C.I.  given  in 
this  section  by  (l,u]  represents  the  formula  used  to  calculate 
interval  estimates. 

1.  C.I.  for  Normal  jJL  with  CT^  Known  [Ref.  2:  pp .  77-80] 

p 

Assuming  that  X  is  distributed  N  ( jJi  >  (J  )  and  using 
(X- jil)  /  (  CT/^n)  as  the  pivotal -quanti  ty ,  a  1007*  C.I.  is 
constructed  thus: 

p[za-y>/2  ■*  <*-/*>/<  cr/./ -=  z(i+'j')/2j  •  y. 

Substituting  c  =  2(i  +  y)/?  *  'z(l-7)/',  and  simPlifying  we  have 

(l,u)  =  (x  -  zCJ/^a  ,  x  +  zCJ/-yTn). 


2 .  C.I.  for  Normal  JJL  with  (J ^  Unknown  [Ref.  2 :  p .  80; 
Ref.  3:  p.  277;  Ref.  1:  p.  381] 

Assuming  X  is  distributed  N and  using 

(X-/I )  /  (s/^Tn)  ,  which  is  distributed  t(n-l),  as  the  pivotal- 

quantity,  a  lOO'T'l  C.I.  is  constructed  thus: 


p  ^u-yvz 

Substituting  t 
fying  we  have: 


(n-l)  -a  (X-  p.)/{s/<yjn)  «:  t(i  +  y ) /2  J  =T‘ 

=  t(l  +  7')/2(n‘1)  =  -^1-7)72 


(n-l)  and  simpli- 


or 


3.  C.I.  for  Bernoulli  p  [Ref.  4:  pp .  376-381] 

The  pivotal -quantity  method  does  not  work  for  the 
Bernoulli  case  and  another  method  will  be  briefly  developed 
here.  This  method  starts  with  two  numbers  a  and  b  such  that 


«  7* 


From  these  limits  we  have: 


P  [X  “  a]  “  (l-'Y)/2  =  P  [x  S  bj  ,  or  equivalently, 

p  [Zxi  S  naj  =  (l-7)/2  =  P  £x.  ss  nb  . 

Now  is  distributed  Binomial  (n,p)  which,  for  k  <  n,  can 

be  explicitly  related  to  the  incomplete  Beta  function  and 
hence  to  the  F  distribution  to  yield: 


u  = 


(n%  *  l)F(1+-y)/2(v1,Y2) _ 

{n  -  nX)  +  (nX  +  1 )P(1+ y )/2^vi » v2) 


nX 

1  - - - - 

(nX  +  (n  -  nX  +  1  ^F(1  +  y  )/2*v2+2’ vi“2^ 

where  v^  =  (2nX  +  2) ,  and  v =  (2n  -  2nX) 

1  and  u  form  a  lOO'/li  conservative  random  interval  thus: 
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p  [i  «  p  « u]  z  y. 


The  outcome  (l,u)  is  a  conservative  C.I.  in  the  sense  that 
for  this  discrete  distribution  the  confidence  that  (l,u) 
contains  p  is  at  least  lOO'y‘3. 

4.  C.I.  for  Normal  (J^  with  /I  Known  [Ref.  3:  p.  275] 
Assuming  that  X  is  distributed  N(  fJL,  o2  )  and  using 
/ CT ) “ ,  which  is  distributed  X^(n)»  as  the  pivotal - 
quantity,  a  100 equal  tails  C.I.  is  constructed  thus: 


X2u-y>/2(n)  -=  Bh -=  Xzu+y )/2 


(n) 


■r 


Substituting  qx  =  X2(  1  -y  )  /  2  (n)  »  and  q2  =  X2(i  +  yj/2Cn) 

and  simplifying  we  have: 


P  [l(Xi-/i)2/q2  *<  <J2  c  I^-jLZ^/qJ  =  T,  or 
( 1 » u )  =  (  Y.xi2  ~  n^2)(l/q2  ,  l/q^. 


5 .  C.I.  for  Normal  _sL  with  fj.  Unknown  [Ref.  1:  p.  382; 
Ref.  3:  p.  277] 

2 

Assuming  that  X  is  distributed  N ( JJL ,  d  )  and  using 

2  2  2 
(n-l)s  /  (J  ,  which  is  distributed  X  (n-1) ,  as  the  pivotal- 

quantity,  a  lOO'Y't  equal  tails  C.I.  is  constructed  thus: 


p[x2u-7)/2 !n-l)  «  (n-l)s2/<72-»  X2(l.y)/2(n-l>]  -  7- 


.. 
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Substituting  qx  =  X  2  (1  ~y  )  /  2  (n"  1')  ’  and  q2  ~  X  (l  +  y)/21'n  3" 
and  simplifying  we  have: 


P  £(n-l)s2/q2  G2  **  (n-l)s2/q1  =  7*. 

(l,u)  =  (n-l)s2(l/q2  ,  l/q^- 


or 


kl! 


5_  c.I.  for  Exponential  X  0r  M  [^ef-  3:  P*  ^79] 
Assuming  that  X1,X2,...XJl  are  exponential  random 
variables  with  parameter  X  and  using  2 AnX,  which  is  dis¬ 
tributed  X2{ln),  as  the  pivotal -quantity ,  a  100T%  equal 
tails  C.I.  is  constructed  thus : 


X2(i-y)/2(2n)/2nX  "e2^nX  "*  X2(i+y)/2(2n)  =^‘ 


Simplifying  we  have 


X2(i-7)/2(2n-)/2nX  X2{i+y)/2(2n)/2nX^ 


or 


( i » u)  =  (X2(i-y)/2(2n)/2nJ  •  X2fi+y)/2(2n)/2nx)- 

The  C.I.  for  the  mean  time  to  failure  (  JJ- =  l/\ )  is 
constructed  by  inverting  the  above  interval  to  yield: 

<l,u)  =  (2nx/X2(i+y  )/2(2n)  ’  2n5E/X2(l-7  )/2(2n)) 

[Ref.  4s  p.  382]  . 
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B.  THEORY  FOR  TWO -POPULATION  CONFIDENCE  INTERVAL  ESTIMATION 
The  confidence  interval  estimates  for  two -population 
situations  are  discussed  here  in  the  same  order  as  they 
appear  in  the  User  Guide  for  Program  2  (Appendix  A) .  All 
of  these  estimates  use  the  pivotal -quantity  method  discussed 
earlier  [Ref.  1:  pp .  379-389].  As  in  the  one-population 
case  above,  the  C.I.  given  by  (l,u)  represents  the  formula 
used  to  calculate  the  interval  estimate. 

1.  C.I.  For  Bernoulli  P^"Py  f°r  Large  m  and  n  [Ref.  2: 
p.  249] 

Large  m  and  n  means  npy,  mpY,n(l-p^),  and  m(l-py)  all 
greater  than  five.  With  this  condition  met  and  assuming  that 
X  and  Y  are  independent  and  normally  distributed, 

(X  -  Y)  -  (px  -  pY) 

Vsx^/ n  +  sY2/m 


which  is  distributed  approximately  N(0,1),  is  used  as  the 
pivotal -quantity .  A  100  7»  C.I.  is  constructed  thus: 


-z 


(i+7V2 


(X  -  Y)  -  (px  -  py) 


-/x(l-X)/n  +  Y(l-Y)/m  (1  +  TV2 


=  7- 


Substituting  c  =  '/xTl  -X)/n  +  Y(l-Y)/m  ,  and  simplifying  yields 


P 


X~  Y  "  cz(l  +  Y  )/2  ■C'PX"PY 


x-y*cZ(i+Y)/21  y. 


or 
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(l,u)  =  (x-y  -  cz^1  +  yy2  ,  x-y  +cz(l+yy2). 


2.  C.l.  for  Normal  /VMy  for  X  and  Y  Paired  [Ref.  2: 
p.  123] 

Assuming  that  X  and  Y  are  normally  distributed  and 
letting  D  =  X  -  7,  and  using  (D-  Md)/cv  ■^Jn)  ,  which  is 
distributed  t(n-l),  as  the  pivotal-quantity,  a  100^^  C.l.  is 
constructed  thus: 


t(l-7')/2(n‘l)  ^  t(i  +  y  )/2^n~1  ^  ~Y' 


Substituting  Jl^  *  [JL^  -  JULY>  and  t  =  t  ^1  +  y  ^  /  2  (n- i ) 
-t  (i-V  )/■?  (n and  simplifying  yields: 


P 


tV  fn  ***  Mx-/^Y 


D  +  tSD/ 


or 


Id^-CEd^Vn 

(l,u)  =  (d  -  tSp/ y  n,  d+tSjyyn)  >  where  sD=  - (n^T5 - 

and  d£  =  Xj,  -  yi. 

3.  C.l.  for  Normal  fly-- fJiY  with  (Jx2  =  (Jy"  =  (J 2 
[Ref.  2:  p.  125] 

Assuming  that  X  and  Y  are  independent  and  using 
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-/ 


,  which  is  distributed  t(m+n-2). 


H  - Y)  -  [Mx  -  Ike \ 

(l/o  ♦  l/n)S  2 

as  the  pivotal -quantity .  Here  and  elsewhere  in  this  paper, 
n  and  m  represent  the  number  of  data  points  in  X  and  Y, 
respectively.  A  100'yi  C.I.  is  constructed  thus: 


(X-Y)  -  (  Uy-Uy) 

'V )  /o(m+n-2)  *=  - 5-  -=  tfl  /„(m+n-2) 

(1"/)/2  (l/m+l/n)Sp2  (l+7)/2 


=  7 


Substituting  t  =  t^ty^,(ni+n-’)  =  -t^_y  ^  ,(m+n-2)  and 
simplifying  yields: 


X  -  Y  -  t  VUA+  l/n)S  2  JHX~  X  -  Y  +  t 


=7 


or  ( 1  * u )  =  (x  -  y  -  t  y(l/m  +  l/n)s  2  ,  x -y  +  t  +  l/n)s  2) , 

.r'  .H 


4. 


where  s 


-  x)2  +  ^(yi  -  y)' 

(m  +  n  -  2 ) 


2  2 

C.I.  for  Normal  fjLx~  JJ.  Y  with  ^X  and  (Jy  Known 


[Ref.  2:  p.  123] 

Assuming  X  and  Y  are  independent  and  using 


(x  -  y)  -  (  n  -n  ) 

- p=  v  ==  g=== —  ,  which  is  distributed  N(0,1),  as 

V  Ox  /n  +  O'y  /m 

the  pivotal -quantity ,  a  10o7^  C.I.  is  constructed  thus: 
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P  *  „  (X  -  Y )  -  (  Uy  -Liy) 

P  "  (l  +  T)/2  x  - -  ■  ■  -  .  <  z 

YCTx2/n  +  0^f2/m 


(i+r)/2  =r- 


Substituting  z  *  z 
have : 


(1  +  7)72  =  'Z(l-7)/2 


and  simplifying  we 


X  -  Y  -  z  V  CTx2/n  +  (7Y2/m  <  /£x-  /Zy<  X  -  Y  +  z  \  (Jx2/n  +  (TY2/m  =7, 


or  (l,u)  =  (5- y-  z  V(Jx2/n+  <Jy2/m  ,  x  -  y  +  z  VcTx2/n  +  cry2/m) 


S.  C.I.  for  Normal  (Ty2/ <Jy2  [Ref.  4:  p.  464] 

Assuming  X  and  Y  are  independent  and  using 
7  2  7  2 

(Sy  /  CTy  )/(SxV  (Jx  ),  which  is  distributed  F(m-l,n-l),  as 
the  pivotal -quantity ,  a  1 0 0 7 7  C.I.  is  constructed  thus: 


3y2/aY2 


P  F(l-r)/2(m'1,n-1)  P(l-T)/2<n'1,n*1) 

X  X 


Substituting  F(1.T)/2(m-l,„-l)  -  .  and 


simplifying  we  have: 


s  2/s  2  tfx2 

pd+7)/2^'1,m"1)<!  ay2  ^  (Sx2/sx2)F(i+7)/2{m_1’n"1)  =  7- 


2/2 


or  (l,u) =  ( 


F(i  +  yy2(n-i,m-1)  ’  (sx2/sy2)F(l+y  ^2(m-l,n-l)  ) 


6.  C.I.  for  Exponential  Xx^Xy  =  A y^M  x  tRef-  4: 


p.  466] 


Assuming  X  and  Y  are  independent  and  using  X  Xx^Xv 
which  is  distributed  F(2n,2m),  as  the  pivotal -quantity ,  a 
100  To  C.I.  is  constructed  thus: 

P  2^2n,  2m)  •<  X  \X/Y  \y  “<  +  ^  ^2(2n,  2m)  =T* 

Simplifying  we  have: 

P  (Y/x)F(l_y  )/2{2n,2m)  -c\x/XY  <  (Y/X)F(l+y  )/2(2n,2m)J  =  J. 

or  (l,u)  =  (  {  Y/X)F(l_yy2(2n,2m)  ,  (Y/x)F(1+y  y2(2n,  2m) ) . 


C.  THEORY  FOR  ONE-POPULATION  HYPOTHESIS  TESTS 

Three  types  of  hypothesis  tests  are  performed  in  Program  3; 
these  are  upper-tailed ,  lower-tailed  and  two-tailed  tests.  The 
general  procedure  used  for  all  three  tests  is  the  same. 
Basically,  a  test  statistic  T,  which  has  an  assumed  distribu¬ 
tion,  is  computed  with  user-supplied  data  and  then  compared 
to  a  critical  region  defined  by  Tq .  This  TQ  value  is  deter¬ 
mined  by  the  type  of  test,  the  assumed  distribution,  and  the 
user-supplied  CC  value.  For  two-tailed  tests  there  are 
actually  two  values  of  Tq  used;  these  values  will  be  denoted 


here  by  T1  and  T 2*  In  performing  the  test  then,  if  the  test 
statistic  T  is  outside  the  critical  region,  then  we  accept 
our  original  hypothesis  Hq  ;  otherwise,  we  reject  Hq  and 
accept  the  alternative  H^. 

The  following  graphs  illustrate  the  three  types  of  tests 
discussed.  These  graphs  represent  a  probability  density 
function,  f(T),  where  the  shaded  area  under  each  curve 
represents  values  from  a  corresponding  cumulative  distribution 
function,  F(T). 


Lower-Tailed  Test 


Two-Tailed  Test  for  Unsymmetric  Distributions  (Chi-Sauare.  F) 
f(T)i 


Median 


An  upper-tailed  test  is  performed  whenever  the  user 
enters  a  +1,  choosing  the  upper-tail  alternate  hypothesis, 
during  step  one  of  the  solution.  A  lower-tailed  test  requires 
the  user  to  enter  a  -1  during  the  solution,  and  a  two-tailed 
test  requires  a  0  entry.  For  two-tailed  unsymmetric  tests 
the  relationship  between  the  test  statistic  and  the  median 
determines  which  tail  is  actually  used  for  the  test.  If 
the  test  statistic  is  greater  than  the  median  then  the  two- 
tailed  test  is  performed  using  the  upper  tail  Tq  value.  For 
test  statistic  values  less  than  the  median  the  test  uses  the 
lower  tail  To  value.  However,  there  is  no  general  agreement 
that  a  two-tailed  test  should  be  performed  this  way. 

The  test  statistics  used  in  Program  5  will  now  be  listed 
in  the  same  order  in  which  they  appear  in  the  User  Guide. 

This  listing  will  also  include  a  reference  where  more  informa¬ 
tion  concerning  a  particular  subject  can  be  found.  Additionally, 
the  assumed  distribution  for  each  test  statistic  will  appear 
with  that  test  statistic. 

1.  Test  Statistic  for  Normal  with  (J  Known 
(Ref.  1:  p.  431] 


,  using  N(0 , 1) 


2.  Test  Statistic  for  Normal  Uo  with  o'2  Unknown 
[Ref.  1:  p.  431] 


30 


v 

t  i 


I 

I 

f  » » 


Li 


T  -  T.  -  X>i  *  -  VlH 

1  -  1  _  9  > 


using  XZ^n) 


For  n  >  65 


T  = 


T '  -  n 


2  n 


using  N(0,1) 


? 

5.  Test  Statistic  for  Normal  (JQ  with  ^  Unknown 


[Ref.  1:  p.  432;  Ref.  2,  p.  104] 
For  n  <  64 


31 


(n-1 )sv2 


using  X  (n-1) 


For  n  >  64 


T'  -  (n-1 ) 
2 (n-1 ) 


using  N(0 , 1) 


6.  Test  Statistic  for  Exponential  JJ.  =  l/^0  [Ref. 
p.  279] 

For  n  <  32 


T  =  T* 


Mo 


using  X2(2n) 


For  n  >  32 


T  = 


using  N(0,1) 


7.  Test  Statistic  for  Poisson  X0  [Ref-  2:  p.  248] 


For  n  <  30 

T  =  n\0  . 


2 


using  X  (2nx) 
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H 


For  n  >  30 


using  N(0 , 1) 


D.  THEORY  FOR  TWO -POPULATION  HYPOTHESIS  TESTS 

The  two-population  hypothesis  tests  are  performed  in  the 
same  manner  as  the  one -population  tests.  For  an  explanation 
of  these  tests,  refer  to  C  above.  The  test  statistics  used 
in  Program  4  will  now  be  listed  in  the  same  order  in  which 
they  appear  in  the  User  Guide.  This  listing  will  include 
the  assumed  distribution  and  applicable  references. 

1.  Test  Statistic  for  Bernoulli  P^  =  Py  [Ref.  2:  p.  249] 


(x  -  y) 


T  = 


-  WWj  (^W)  5 +  i) 


,  using  N (0 , 1) 


2.  Test  Statistic  for  Normal  jj, ^  =  JJLy  for  X,Y  Paired 


[Ref.  2:  p.  121] 

T  =  cTC/n/s^)  ,  using  t(n-l)  for  n  <  31,  and 

N  (0 , 1)  for  n  >  31 
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3.  Test  Statistic  for  Normal  for  X  and  Y 


Independent  with  =  (Jy“  1:  PP*  434-435; 

Ref.  2:  p.  116 


,/T+  I 

p  y  n  m 


,  where  s  = 
P 


(\n-l  )sx2  +  (m-1  )sy2 
n  +  m  -  2 


and  using  tCn+m-2)  for  n  <  32,  and 
N(0 , 1)  for  n  ^  32 

4.  Test  Statistic  for  Normal  jj_  ^  =  ^Y  for  X  and  Y 


Independent  with  ,  (Jy  Known  [Ref.  2:  p.  119] 


_ x  -  y 

T  =  ( - 2 - 2  ’  using 

“\  £jl.  +  _5x_ 

V  n  m 


5.  Test  Statistic  for  Normal  Mx  =  My  for  X  and  Y 


Independent  with  (Jy  f  (JY  [Ref.  2:  p.  119] 


L  2  _  2 

5X  +  Y 
n  m 


,  using  N(0,1)  for  df  >_  30,  and 
t (df )  for  df  <  30,  where 


Independent  [Ref.  2:  p.  Ill] 


T  =  — — ,  using  F(n-l,m-l) 


7.  Test  Statistic  for  Normal  P  =  0  [Ref.  1:  pp .  492- 
493;  Ref.  2:  p.  200] 


using  t(n-2)  for  n  <  28,  and 
N  CO  ,  1)  for  n  >  28 


8.  Test  Statistic  for  Exponential  =  Xy  ^or  ^  an^  ^ 

Independent  [Ref.  4:  p.  310] 

T  =  ,  using  F(2n,2m) 

y 


E.  THEORY  FOR  METHODS  OF  APPROXIMATION 

The  methods  used  to  approximate  distribution  values  for 
this  package  of  programs  come  from  many  sources,  and  any 
complete  discussion  of  the  theory  involved  would  be  beyond 
the  scope  of  this  paper.  Therefore,  while  every  method  which 
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has  been  used  will  be  referenced,  only  those  methods  unique 
to  this  programming  effort  will  be  discussed  in  theoretical 
detail.  With  this  in  mind,  a  combination  listing  and  dis¬ 
cussion  will  follow  which  will  trace  the  approximation  methods 
used  in  this  paper.  The  approximations  unique  to  this  effort 
are  the  inverse  cumulative  distribution  functions  (CDF)  for 
the  Chi-Square,  Student's  t  and  F  distributions  as  well  as 
the  Multinomial  approximation.  The  inverse  Normal  CDF 
approximation  discussed  here  was  used  in  Professor  Zhena's 
programs.  All  other  approximations  used  in  this  paper  are 
modifications  of  methods  from  the  TI-59  Statistics  Module. 

For  a  better  discussion  of  the  exact  methods  used  in  these 
cases,  the  TI-59  Applied  Statistics  Manual  should  be  consulted 
[Ref.  7]. 

There  are  two  types  of  inverse  CDF  approximations  used  in 
this  paper.  The  first  two  programs  use  a  less  accurate 
Type  I  approximation,  while  the  distribution  programs  use 
a  closely  related,  but  more  accurate,  Type  II  approximation. 
What  follows  is  a  discussion  of  the  methods  used  for 
approximating  the  inverse  Normal,  Chi-Square,  Student's  t  and 
F  distrubitons ,  and  the  Multinomial  distribution. 

1.  Inverse  Normal  CDF  Approximation  [Ref.  5:  p.  933] 
a.  Type  I  Approximation 

This  approximation  is  used  in  Programs  1  and  2 
for  inverse  Normal  CDF  values  and  as  a  subroutine  for  the  Chi- 
Square,  t  and  F  approximations  in  those  same  programs.  The 
Type  I  approximation  uses  the  following  set  of  equations  and 


constants  to  approximate  inverse  Normal  CDF  values,  given  the 
input  probability  p.  Programs  1  and  2  have  the  limitation 
that  p  be  greater  than  .5. 


Z 

P 


=  t  - 


c0  +  Clt  +  c2t2 
1  +  dtt  +  d2t2  +  d3t3 


,  where  t 


Vlnd/d-p)2) 


and 


c  0  =  2.S15517 
c1  =  .802855 
Co  =  .010328 

L, 


dx  =  1.432788 
d2  =  .189269 
d3  »  .001308 


b.  Type  II  Approximation 

The  Type  II  approximation  is  almost  identical 

to  the.  Type  I  approximation.  The  only  difference  is  the 

addition  of  a  function  which  removes  the  limitation  that  p 

be  greater  than  .5.  This  function  uses  the  symmetric 

quality  of  the  Normal  distribution  and  returns  the  negative 

of  the  approximation  for  1-p,  whenever  p  is  less  than  .5. 

2.  Inverse  Chi-Square  CDF  Approximation  [Ref.  5:  p.  941] 

a.  Type  I  Approximation 

This  approximation  is  used  in  Program  1  only 

and  uses  the  inverse  Normal  CDF  approximation  described 

above.  The  inverse  Chi-Square  CDF  approximation  uses  the 

2 

following  set  of  equations  and  calculates  both  X 
and  X  (p) (v)  given  the  input  probability  p.  The  reference 
listed  above  limits  degrees  of  freedom,  v,  for  this  approxima¬ 
tion  to  values  above  30;  however,  as  can  be  seen  in  Appendix  B 
values  well  below  v  =  30  produce  acceptable  approximations. 
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X2pO)  =  V(1  -  a  +  Zp  a)J  ,  where  a  =  ^ 

b.  Type  II  Approximation 

This  approximation  is  used  only  in  Program  7  and 
incorporates  the  Type  I  inverse  Chi-Square  CDF  approximation 
with  an  Accuracy  Enhancing  Technique  (AET)  which  requires  the 
highly  accurate  forward  Chi-Square  CDF  approximation  contained 
in  that  same  program.  This  AET  involves  taking  the  output  of 
the  Type  I  approximation  and  using  it  in  the  forward  CDF 
approximation  to  obtain  an  estimate,  £ ,  of  the  actual  proba¬ 
bility  achieved  by  the  Type  I  approximation.  This  estimate 
is  then  used  to  correct  the  inverse  approximation  input,  p, 
for  any  difference  between  desired  and  actual  probability. 

The  corrected  inverse  input,  p',  is  computed  using  the  follow¬ 
ing  formula: 

P'  ■  P  +  (P  -  P)  • 

By  using  p'  as  the  new  input  for  the  inverse  CDF  approxima¬ 
tion,  a  more  accurate  approximation  is  achieved  (see  Appendix  B) 
This  AET  is  not  used  for  v  >  30  since  the  Type  I  approximation 
is  then  quite  accurate. 

3.  Inverse  Student's  t  CDF  Approximation 
a.  Type  I  Approximation 

This  approximation  is  used  in  Programs  1,  2,  8 
and  9  and  requires  the  inverse  Normal  CDF  approximation 
already  described.  The  following  equation  represents  one  of 
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several  approximations  to  the  inverse  t  developed  by  Professor 
Donald  P.  Gaver ,  Naval  Postgraduate  School.  This  approxima¬ 
tion  is  limited  to  values  of  v  >  2  [Ref.  8] . 


tp(v)  =  Zp  |  1  +  Zp 


2  1-1  + 


10 


3(v  -  1-57) 


v2:2 


b.  Type  II  Approximation 

This  approximation  is  used  only  in  Program  8 
and,  like  the  Type  II  Chi-Square  approximation,  uses  an  AET 
to  increase  the  accuracy  of  the  Type  I  approximation.  The 
procedures  for  this  AET  are  exactly  the  same  as  in  the  Chi- 
Square  approximation.  Additionally,  the  restriction  on  v 
has  been  removed  by  using  the  following  relationship  to 
generate  inverse  approximations  when  v  =  1. 


t„(l)  = 


tan(l-p) 


,  where  for  p  <  .5,  2p  is  used  in  place  of  p 


4.  Inverse  F  CDF  Approximation  [Ref.  5:  p.  947] 
a.  Type  I  Approximation 

This  approximation  is  used  in  Programs  1  and  2 
and  uses  the  inverse  Normal  approximation  discussed  earlier. 
The  following  set  of  equations  is  used  to  generate  the 
approximation : 


y-i 


where  w 


Z-p  ^(h  +  k) 
h 


h 


*  2(i/(v1  -  i )1  +  "i/(v2  -  ITT  ’ 


and  v^/l  *  v2^1 


b.  Type  II  Approximation 

This  approximation  uses  the  AET  discussed  earlier 
along  with  two  more  techniques  to  provide  inverse  approxima¬ 
tions  for  all  values  of  v^  and  v^ .  For  the  case  where  v^  or 
v7  -  1 ,  the  Type  I  inverse  t  approximation  discussed  earlier 
is  used  in  the  following  relationships  to  generate  the  required 
inverse  F  CDF  approximations: 


F  (l,v)  =  (t(l+p)/2(v))2  ,  or  equivalently 


vv'i)  =  • 


For  the  case  where  v^  =  ~  1,  the  following  relationship 

is  used: 

Fp^  ,l  *  =  tan(l-p)  • 
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5 .  Multinomial  Approximation 

The  method  used  to  generate  multinomial  density  values 
is  not  truly  unique;  however,  it  uses  the  following  equation 
in  a  way  that  minimizes  rounding  errors: 


( nl ’ 


n 


2’ 


,nk)  = 


N 


(n^ ! 


)(n2!  ). .  .(nv!  )(P1nl)(P2n2^- •  -{Pknk) 


Computations  are  accomplished  in  the  following  order  to 
avoid,  as  much  as  possible,  multiplying  extremely  small 
values  by  extremely  large  values  : 


4(Plnl)\  /{N  -  l)(p2n2)N 


V 


■2' 


(N  -  k)(pknk)'' 


n. 


k* 


(N  -  k  -  1)! 
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IV.  TI-59 


The  choice  of  the  TI-59  as  the  calculator  for  this  pro¬ 
gramming  effort  was  based  on  two  factors.  Professor  Zehna's 
original  programs  were  written  for  the  TI-59;  and  each  student 
in  the  Operations  Research  (OR)  curriculum  is  issued  a  TI-59 
for  use  in  basic  probability  and  statistics  courses.  In 
general,  the  use  of  hand-held  programmable  calculators  has 
been  shown  [Ref.  6:  p.  1]  to  increase  student  learning  and 
capability.  Further,  it  is  intended  that  the  programs 
described  in  this  paper  will  be  used  by  OR  students  in  their 
coursework . 

Using  the  TI-59  offered  some  disadvantages  and  some 
advantages  in  developing  the  programming  package  presented 
here.  The  disadvantages  of  limited  storage,  slo\^  computation 
time  and  awkward  data  entry  sequence  are  discussed  elsewhere 
in  this  paper.  There  are  two  advantages  of  the  TI-59,  how¬ 
ever,  that  deserve  noting  here.  The  programming  steps  and 
procedures  used  in  the  TI-59  are  easily  learned  and  logical. 
This  ease  of  programming  makes  complex  computational  methods 
easy  to  program.  Another  advantage  of  the  TI-59  is  its 
compatibility  with  the  PC-100C  printer.  Using  the  printer/ 
calculator  combination  greatly  simplifies  writing,  editing, 
and  error  diagnosis.  These  advantages  of  the  TI-59  make 
programming  relatively  easy  and  also  allow  the  capabilities 
of  the  TI-59  to  be  used  more  fully. 
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V.  ALTERNATE  SOLUTIONS 


The  problems  presented  in  this  paper  could  have  been 
solved  in  any  number  of  equally  valid  ways.  This  section 
will  briefly  discuss  the  specific  alternatives  which  could 
increase  the  capability  of  the  solutions  employed  here. 

This  discussion  will  first  focus  on  alternative  hardware 
which  could  be  used,  and  then  on  program  changes  which  might 
be  made. 

A.  HARDWARE  ALTERNATIVES 

As  discussed  earlier,  the  use  of  the  TI-59  calculator  is 
appropriate  for  this  programming  package;  however,  the  methods 
and  techniques  used  in  this  paper  are  equally  suited  to  other 
programmable  calculators  or  computers.  Indeed,  the  use  of 
any  other  calculator  or  computer  with  more  storage  capability 
than  the  TI-59  might  be  a  better  vehicle  for  the  package  of 
programs  presented  here.  These  programs  require  a  total  of 
14  magnetic  cards  (22  separate  sides)  using  the  TI-59. 

Current  console  model  computers  could  easily  store  all  of 
these  programs  at  the  same  time. 

The  TI-59  has  the  ability  to  use  program  modules,  and 
the  nine  programs  from  this  package  could  easily  be  combined 
to  form  the  basis  of  a  new  module.  By  careful  elimination  of 
redundant  functions,  such  a  module  could  also  accommodate  an 
ANOVA  and/or  regression  package.  A  module  of  this  type  might 
find  widespread  acceptance,  especially  in  the  classroom. 
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B.  PROGRAMMING  ALTERNATIVES 


The  choice  of  distribution  approximating  methods  used 
here  and  the  general  layout  of  the  entire  package  might 
both  be  improved  with  some  additional  effort.  The  methods 
of  approximation  used  here  were  chosen  rather  arbitrarily. 

It  is  possible  that  more  appropriate  methods  of  approximation 
exist.  More  appropriate  methods  might  include  methods  with 
fewer  program  steps  and  equal  accuracy,  as  well  as  methods 
which  take  less  computation  time  with  equal  accuracy. 

The  current  layout  of  the  nine  programs  in  this  package 
might  be  improved  by  reducing  the  number  of  programs.  The 
present  program  size  precludes  this;  however,  by  eliminating 
some  program  functions,  it  might  be  possible  to  organize  all 
of  the  inverse  CDF  approximations  in  one  program.  Shorter 
approximating  methods  might  make  possible  a  similar  program, 
but  with  both  inverse  CDF  and  regular  CDF  approximations. 

The  extensive  user  guides  with  their  sample  problems  are 
useful  for  students,  but  could  prove  awkward  for  a  more 
experienced  user.  An  obvious  alternative  would  be  a  shortened 
user  guide  directed  at  those  more  familiar  with  the  programs. 


VI.  CONCLUSION 


Two  problems  were  presented  for  solution  in  this  paper. 

One  problem  involved  expanding  and  modifying  an  existing  set 
of  TI-59  programs  into  a  user-friendly  package.  A  second 
problem  involved  developing  a  set  of  distribution  approximating 
programs.  The  solution  to  the  first  problem  incorporated 
increased  capability,  standardized  data  entry,  and  detailed 
user  guides  into  a  package  of  nine  programs.  The  first 
priority  in  this  solution  was  providing  a  format  compatible 
with  a  student's  needs  while  maintaining  the  capability 
required  by  a  more  experienced  user.  The  TI-59  calculator 
proved  to  be  a  very  useful  tool  in  this  solution,  and 
demonstrated  the  generally  unused  capability  of  the  current 
generation  of  hand-held  programmable  calculators. 

The  approximation  programs  presented  as  a  solution  to  the 
Second  problem  mentioned  above  provide  accurate  and  compre¬ 
hensive  approximations.  These  programs  practically  eliminate 
the  need  for  tables  of  values  and  solve  the  interpolation 
problem  present  in  all  such  tables. 

Together  these  two  solutions  represent  a  package  with 
considerable  capability  in  computing  confidence  intervals, 
performing  hypothesis  tests,  and  generating  approximate  distri¬ 
bution  values.  A  comprehensive  set  of  user  guides  makes 
this  same  capability  available  even  to  inexperienced  users. 

The  methods  used  in  preparing  this  TI-59  programming 
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package  are  directly  applicable  to  other  calculators  or 
computers . 


APPENDIX  A 


PROGRAM  1  USER  GUIDE  -  One -Population  Confidence  Intervals 
INTRODUCTION;  The  purpose  of  this  program  is  to  compute  100^% 
confidence  intervals  (l,u)  or  bounds  £l  and  uj  for  the  following  one- 
population  situations; 

NORMAL  fj.  with  (J 2  known 

NORMAL  fj.  with  (J^  unknown 

BERNOULLI  p 

NORMAL  (J2  with  jj.  known 

NORMAL  CT2  with  jj.  unknown 

EXPONENTIAL  \ot  j±  . 

The  routines  in  this  program  require  percentiles  from  either  the 
Normal,  Chi-Square,  Student's  t  or  F  distributions.  Each  routine  will 
automatically  generate  an  approximate  percentile;  however,  when 
additional  accuracy  is  desired  or  small  sample  sizes  are  involved  the 
use  of  percentile  values  from  either  standard  tables  or  the  distribution 
approximating  programs  is  recommended.  In  step  two  of  each  routine  the 
user  can  choose  to  accept  the  approximate  percentile,  by  storing  the 
appropriate  percentage  in  storage  register  09  (Rq^) <  or  he  can  store 
the  percentile  value  in  R^  .  Some  routines  also  require  percentile 
values  in  R^ . 

GENERAL  PROCEDURES; 

1.  Use  any  library  module,  and  after  reading  all  three  card  sides, 
press  |  D  '  to  repartition  (639-39)  • 

2 .  For  data  entry  press  m  followed  by  data  point  x^,  [r/s| ,  x^, 
R/S  .  etc.  for  each  x^  (i=l ,2, . . .  ,n)  until  all  points  have  been 
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entered.  Mistakes  in  data  entry  should  be  corrected  immediately  by  reen¬ 
tering  the  unwanted  point  and  pressing  INV  2nd  Xt  ,  then  enter  the 


correct  data  point  and  press 


.  Alternate  data  entry  using 


summary  statistics  is  detailed  in  applicable  routines. 

3.  For  one-sided  confidence  bounds  rather  than  intervals  replace 
(1  +y)/2  with  y  and  (l-y)/2  with  l-*/  everywhere  they  appear 

(e.g.  z(l+y)/2  becomes  zyS>  and-  Proceed-  as  usual,  ignoring  1  or  u 

as  appropriate. 

4 .  When  solving  consecutive  problems ,  care  should  be  taken  to  clear 
all  previously  used  registers.  Pressing  m  will  clear  all  registers . 
PROGRAM  1  SPECIFIC  PROCEDURES: 


C.I.  For  NORMAL  J±  with  (J2  known 

1.  Enter  data  using  [  D  |  and  store  (J  in  RQ^  (Alternate  entry: 
store  n  in  RQ^,  x  in  Rgg  and  (J  in  Rq^) 

2.  Store  either  (l+')'')/2  in  OR  z^i+y)/2  in  R11 

3 .  Press  m  1  is  displayed,  then 
press  f  xgt (  u  is  displayed 

C.I.  For  NORMAL  ^  with  (J2  unknown 

1.  Enter  data  using  j  D  |(Alternate  entry:  store  n  in  r03' 
s  in  RQr,  and  x  in  Rqq) 

2.  Store  either  ( l+y)/2  in  Rq^  OR  b(i+'^')/2^n_^'^  in  R11 

3.  Press  j  C  |  1  is  displayed,  then 
press  fx£t |  u  is  displayed 


.  V 


* 


G.I.  For  BERNOULLI  P 

1 .  Store  nx  in  Rq^  and  n  in  Rq^ 

2.  Store  either  (l+^)/2  in  Rq^  OR  ^y'2(2nx+2f2n-2nx)  in  R^ 


AND 


F^1+<y)/2(2n-2nx+2,2nx)  in  R^ 

3 .  Press  m  1  is  displayed,  then 
press  |x%t [  u  is  displayed 
G.I.  For  NORMAL  ( J 2  with  known 

1 .  Enter  data  using  [  D  |  and  store  j_l  in  Rqq  (Alternate  entry : 

store  ^x^2  in  Rq^,  n  in  Rq3  anci  M  in  RC8^ 

2.  Store  either  (1+7')/2  in  R^  OR  X^l-y)^11)  in  R11 


AND 
,2 


3 .  Press  mn  1  is  displayed,  then 
press  (x%± |  u  is  displayed 
G.I.  For  NORMAL  2  with  unknown 


X  (i+7)/2(n)  in  R13 


1.  Enter  data  using  |  D  j  (Alternate  entry :  store  n  in  R^  and 
(n-l)s2  in  R12) 

2.  Store  either  ( l+y)/Z  in  RQ^  OR  X2(l-'y)/2^n~^^  in  R11 

AND 

X2(l+7)/2(n'1)  ln  R13 

3«  Press  HD  1  is  displayed,  then 
press  ]  xVt  |  u  is  displayed 
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G.I.  For  EXPONENTIAL  \  or  fj* 

1 .  Enter  data  using  cn  (Alternate  entry t  store  n  in  and 
x  in  Rq8) 

2.  Store  either  (1  +'Y')/2  in  R^  OR  ^_~yy^(2n)  in  R^ 

AND 

X2(l+7)/2<2”>  ln  S13 

3 .  Press  HD  1  is  displayed,  then 
press  [x£t |  u  is  displayed 

*  4 .  To  compute  a  confidence  interval  or  bound  for  V\  press 

QD  rather  than  m  above . 

PROGRAM  1  ADDITIONAL  CAPABILITIES! 


Inverse  Normal  CDF  Approximation 


1 .  Store  p  in  R^  (p 


•5) 


2 .  Press  z  is  displayed,  and 


P 

-J.O  ®  VCLJ._LCJ.U_LC  XU  ^ 

Inverse  Chi-Square  CDF  Approximation 

1.  Store  p  in  RQ^  (p  >  .5),  and  store  v  in  R^  (v  /  1) 

2.  Press  JsBR  |  |  X2  |  X2p(v)  is  displayed,  and 

/2 


X^(i-p)(v)  is  available  in  R^ 


1. 

Store 

p  in 

R09 

(p  =»  .5),  in  R20  and  v£  in  R£1 

2. 

Press 

SBR 

RCL 

!  F^(v^ ,v£)  is  displayed 

Inverse 

Student ’ s 

t 

CDF  Approximation 

1. 

Store 

p  in 

R09 

(P  =*  -5)  and  v  in  R^  (v  ^  1) 

2. 

Press 

SBR 

! 

STO 

[  t  (v)  is  displayed 
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PROGRAM  1  LABELS  USED: 


A 

A' 

SBR 

RCL 

3 

3' 

GTO 

STO 

G 

C' 

x*t 

D 

D' 

X2 

E 

E' 

= 

PROGRAM 

1  STORAGE  REGISTER  i 

CONTENTS 

s 

00 

clear 

15 

used 

01 

& 

16 

clear 

02 

l\z 

17 

clear 

03 

n 

18 

clear 

04 

used 

19 

V 

05 

clear 

20 

V1 

06 

clear 

21 

v2 

07 

s  or  (J 

22 

clear 

08 

x  or  fl 

23 

clear 

09 

(l+7)/2  or  ? 

24 

used 

10 

used 

25 

clear 

11 

CDF  value 

26 

clear 

12 

(n-l)s2 

27 

used 

13 

CDF  value 

28 

used 

14 

V 

29-49  clear 

f 


PROGRAM  1  SAMPLE  PROBLEMS; 

1.  Suppose  x  =  69. 7 i  n  =  8,  CT^  =  3>5‘  Find  a  9 0 %  G.I.  for 
estimating  the  mean. 

SOLUTION) 

(1)  Store  n=8  in  Rq^>  x  =  &9-7  in  Rqq  and  (J  =  J  3*5  in  Rq9 

(2)  Store  (l+7)/2  =  -95  in  RQ9  OR  s>g5  =  1.645  in  Rn 

(3)  Press  CD  then  [xS,t  j  to  display 

1  =  68.61179492  1  =  68.61193477 

OR 

u  =  70.78820508  u  =  70.78806523 

2.  Given  the  following  observations  from  a  population  with  a  known 
variance  of  15,  find  a  lower  99%  confidence  bound  on  the  mean. 

165  178  160  199  167  145  1  37  1  82  1  92  1  65 

SOLUTION) 

.  j~  '  • 

(1)  Enter  data  using  j_Djseq.uer.ce  and  store  (J  =  -^  1 5  in  RQy 

(2)  Store  7=  .99  in  RQ9  OR  z_99  =  2.326  in 

(3)  Press  then  fx£tj  to  display 

1  =  168.1502816  1  =  168.1512434 

OR 

ignore  u  ignore  u 

3.  Find  a  90%  G.I.  for  the  mean  of  the  following  test  scores. 

35  34  46  20  38  39  32  49  41  25  18  43  51  38  42  29  59  53  27  33 

SOLUTION; 

(1)  Enter  data  using  ca  sequence 

(2)  Store  (l+y)/2  =  .95  in  RQg  OR  ^(19)  =  1-729  in  Rn 

(3)  Press  ED  then  M  to  display 

1  =  33.39962583  1  =  33.39390952 

OR 

u  =  41.80037417  u  =  41.80609048 
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4 .  Find  a  95?S  upper  confidence  bound  for  jj  given  the  following 
informations  n  =  55.  s  =  15  and  x  =  85. 


SOLUTIONS 

(1) 

Store  n  =  55  in  R  ,  s  =  15  in  RQ^  and  x  =  85  in  RQg 

(2) 

Store  y  =  .95  in  RQ9 

OR  t>95(54)  =  1.673  in  Rn 

(3) 

Press  |  G  |  then  j 

x£t  j 

^  to  display 

ignore  1 

OR 

ignore  1 

u  =  88.38398464 

u  =  88.38380911 

5.  Of  1000  people  treated  with  a  certain  drug  200  showed  a  reaction. 
Find  a  90/2  C.I.  for  the  proportion  of  the  sample  population  that 
will  show  a  reaction. 


SOLUTIONS 


(1) 

(2) 


(3) 


Press  m  then  store  nx  =  200  in  RQ1  and  n  =  1000  in  R 
Store  (l+y)/2  =  .95  in  RQ9  OR  F>95(402,1600)  =1.13  in  R^ 

AND 

F>95(1602,400)  =  1.14  in  R 


Press 


then 


to  display 


1  =  .1793622306 

OR 

u  =  .2219574577 


1  =  .1796719191 
u  =  .2211307235 
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o.  Find  a  95$  upper  confidence  bound  for  the  proportion  of  the  sample 
population  if  nine  of  24  treated  were  affected . 

SOLUTION: 

(1)  Press  m  then  store  nx  =  9  in  Rq^  and  n  =  24  in 

(2)  Store  J  =  -95  in  RQ9  OR  F_95(20,30)  =  1.93  in 

AND 

F<95(32,18)  =  2.13  in  R^ 

(3)  Press  CE  then  jxSt  [  to  display 

ignore  1  ignore  1 

OR 

u  =  .5629409878  u  =  .5626822157 

7.  Find  a  95$  G.I.  for  (J^  and  for  (J  given  the  following  observations 
from  a  sample  population  with  fJi=  65. 0  . 

ICO  15  73  46  65  98  79  38  68  85 
SOLUTION: 

(1)  Enter  data  using  f  D  |  and  store  j[_l=  65  in  Rq8 

(2)  Store  (l+y)/2  =  -975  in  RQ9  OR  X2.025^10)  =  3'^3  ln  fill 

AND 

X2.975^10)  =  20  <48  in  R13 

(3)  Press  m  then  [xSt  |  to  display  (for  (J2) 

1  =  428.7608301  1  =  428.8574219 

OR 

u  =  2726.124604  u  =  2702.461538 

Taking  square  roots  the  following  limits  for  (J  are  displayed 

1  =  20.70654076  1  =  20.70887302 

OR 

u  =  52.21230319  u  =  51.985205 
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8.  Suppose  n  =  15.  =  88476  and  f±  =  30-5  •  Find  a  9?.^  upper 

confidence  bound  on  the  standard  deviation . 

SOLUTION » 

(1)  Press  m  then  store  =  88476  in  R^,  n  =  15  in  RQ^ 

and  (1=  30.5  in  RQg 

(2)  Store  y  =  .975  in  RQ^  OR  X2.025^1'^  =  6,26  in  R11 

AND 

X2.975(15)  =  27.49  in  R13 

(3)  Press  □D  then  |x£t  j  and  take  square  root  to  display 

ignore  1  ignore  1 

OR 

u  =  109.2668647  u  =  109.1078035 

9.  Find  a  95%  C.I.  for  (J  given  the  following  observations. 

100  15  73  46  65  98  79  38  68  85 

SOLUTION  x 

(1)  Enter  data  using  CD 

(2)  Store  {Uj)/Z  =  -975  in  RQ9  OR  X2.025^9)  =  2,70  in  R11 

AND 

^2.975^  =  19,02  in  r13 

(3)  Press  m  then  [S  and  take  square  roots  to  display 

1  =  18.54742  1  =  I8.54896609 

OR 

u  =  49.471926 97  u  =  49.23150151 
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10.  Suppose  n  =  45  and  s2  =  %  find  a  9 0 %  lower  confidence  bound  for 

or2. 

SOLUTION t 

(1)  Press  m  then  store  n  =  45  in  Rq^  and  (n-l)s^  =  1584 
in  R12 

(2)  Store  J  =  . 90  in  RQ9  OR  X*  .10^  =  32,5  in  E11 

AND 

X2.90^)  =  5614  in  R13 

(3)  Press  HD  to  display 

1  =  28.10395538  1  =  28.08510638 

OR 

ignore  u  ignore  u 


11 .  Given  the  time  to  failure  of  an  electron  tube  is  an  exponential 
random  variable,  and  the  sum  of  25  times  to  failure  is  25242. 
Find  a  95%  C.I.  f or  \  . 


SOLUTION i 


(1)  Press  |  D  ~|  then  store  n  =  25  in  RQ^  and  x  =  25242/25 
in  Rq8 

(2)  Store  (Uj)/Z  =  .975  in  RQ9  OR  X2.02 5(5°)=31.92  in  Ru 

AND 


(3) 


Press 


m 


X2.9„(50)=  In 


then 


x*.t 


to  display 


1  =  .0006407042 
u  =  .001414878 


OR 


1  =  .000632 2795 
u  =  .0014048015 
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Six  expensive  pieces  of  equipment  had  the  following  times  to 
failure.  Find  a  95%  G.I.  for  the  mean  time  to  failure. 

233-6  3119-0  258-3  1402.7  612.9  2211.2 


SOLUTION ; 

(1)  Enter  data  using  m 

(2)  Store  (l+7)/2  =  .975  in  RQ9  OR  .025(12)  =  4-40  in  % 

AND 

.975(12)  =  23-34  in  R13 

(3)  Press  E  then  [x£t  j  to  display 

1  =  671.6219738  1  =  671.6109683 

OR 

u  =  3578.128804  u  =  3562.590909 


PROGRAM  2  USER  GUIDE  -  Two-Population  Confidence  Intervals: 
INTRODUCTION :  The  purpose  of  this  program  is  to  compute  100 'Y% 
confidence  intervals  (l,u)  or  bounds  and  uj  for  the  following  two- 
population  situations . 

BERNOULLI  p^-p^  ^or  sample  sizes 

NORMAL  y  ■for  *  and.  Y  paired  with  (J^  unknown 

NORMAL  -or  X  and  Y  independent  with  (Jy2=(J2 

NORMAL  Y  *or  ^  X  Indep  with  and  known 

NORMAL  CJx2/CTy2  dor  X  and  Y  independent 

EXPONENTIAL  X^/\Y  =  MY^Mx  ^or  X  and  Y  independent 

The  routines  in  this  program  require  percentiles  from  either  the 
Normal,  Student's  t  or  F  distributions.  Each  routine  will  automat¬ 
ically  generate  an  approximate  percentile;  however,  when  additional 
accuracy  is  desired  or  small  sample  sizes  are  involved  the  use  of 
percentile  values  from  either  standard  tables  or  the  distribution 
approximating  programs  is  recommended.  In  step  two  of  each  routine  the 
user  can  choose  to  accept  the  approximate  percentile,  by  storing  the 
appropriate  percentage  In  Rq^i  or  he  can  store  the  the  percentile  value 
in  Rjj .  Some  routines  also  require  percentile  values  in  R^. 

GENERAL  PROCEDURES i 

1.  Use  any  library  module,  and  after  reading  all  three  card  sides, 
press  m  to  repartition  (719.29). 

2.  Data  entry  for  Independent  data  (DEI  Sequence),  press  m 
followed  by  data  point  x^.,  r/s  ,  xi+1 ,  R/s  ,  etc.  for  each  x±  (i=l,2, 
...,n).  When  all  of  the  x^'s  have  been  entered  press  HD  followed  by 
data  point  y^  r/s  ,  yi+1 ,  r/s  ,  etc.  for  each  yi  (i=l ,2, . . . ,m) . 
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When  all  of  the  y^’s  have  been  entered  press  [SBR |  GTO  .  Mistakes  in 
data  entry  should  be  corrected  immediately  by  reentering  the  unwanted 
data  point  and  pressing  INV  [2nd}  [13  ,  then  enter  the  correct  data 
point  and  press  R/S  R/S  •  Alternate  data  entry  using  summary 
statistics  is  detailed  in  applicable  routines . 

3.  Data  entry  for  paired  data  (DEP  Sequence),  press [  D  | followed 
by  x^  fxfetj ,  y  ,  [r/s] ,  xi+1>  fxfet  | ,  y±+1 ,  [r/s| ,  etc.  for  each  data 

pair  x^y^  (i=l  ,2 . n) ,  then  press  SBR  (rST|  ■  Mistakes  and 

alternate  entry  are  as  above . 

4.  For  confidence  bounds  rather  than  intervals,  replace  (1 +'y)/2 

with  'Y  everywhere  it  appears  (e.g.  becomes  i'y'(v)  )  aad. 

proceed  as  usual  ignoring  either  1  or  u  as  appropriate . 

5.  When  solving  consecutive  problems,  care  should  be  taken  to 
clear  all  previously  used  storage  registers.  Pressing  m  will  clear 
all  registers. 

PROGRAM  2  SPECIFIC  PROCEDURES; 

* 

C.I.  For  BERNOULLI  p-p,,  for  large  sample  sizes 

— -  '  A  X 

1.  Press  f  D  |  then  store  nx  in  RQ^,  my  in  RQ2  n  in  R m  in  RQ^ 

2.  Store  either  (l+')/)/2  in  R^  OR  z(\+'y)/2  in  R11 

3 .  Press  m  1  is  displayed,  then 
press  [x<t  |  u  is  displayed 

*  Large  here  means  np  ,  mp  ,  n(l-p  ),  m(l-p  )  all  greater  than  five. 

A  X  A  X 


C.I.  For  NORMAL  [lyT  [Ix  for  galr&d'X  and  Y  (n  pairs >  •  *  ■ 

1.  Enter  data  using  DEP  Sequence,  degrees  of  freedom,  v  =  n-1,  will 
be  displayed.  (Alternate  entry;  store  x-y  in  R^Q,  vV*  in  R12 
and  n-1  in  R^) 

2.  Store  either  (1+Y)/Z  in  RQy  OR  t(i+y  ^C11-1)  1x1  R11 

3.  Press  m  1  is  displayed,  then 
press  [x$.t )  u  is  displayed 

4.  When  differences,  x^-y^,  are  given  use  Program  1. 

C.I.  For  NORMAL  with  (j2  =  q^2  =  (J 2  unknown 


1.  Enter  data  using  DEI  Sequence,  degrees  of  freedom,  v  =  n+m-2  is 
displayed.  (Alternate  entry  #1*.  Store  n  in  R^,  m  in  RQ^,  x-y  in 

R10’  Bxi‘x)  ln  ri6*  Edi-r)  in  r26  and  (n+m“2)  in  ri9) 

(Alternate  entry  #2:  Store  n  in  RQ„,  £x±  in  RQ1 ,  ^x±2  in  RQ2, 
then  press  S  ,  store  a  in  RQ^,  £yt  in  RQ1  and  Yj*  ^  111  Ro2’ 
then  press  js3R )  [~GTo] ,  v  is  displayed) 

2.  Store  either  (l+y)/2  in  RQ9  OR  t(i+'y)/2^n+m~2^  in  Rn 

3 .  Press  m  1  is  displayed,  then 
press  [x£t |  u  is  displayed 


C.I.  For  NORMAL  wlth  (J2  and  known 


Enter  data  using  DEI  Sequence,  then  store  (Jl2/n  in  R.  and 
(JY2/m  in  (Alternate  entry:  store  x-y  in  Rin  and 


's/iO/M  +  ( CTY2/m)  in  R12) 

2.  Store  either  (l+y)/2  in  RQ^  OR 
3  •  Press  m  1  is  displayed,  then 
press  [x^t|  u  is  displayed 


09  za.  “(i+y)/2 


10 


in  R, 


11 


&CL 


C.I.  For  NORMAL  (J*/  g-,2 

1.  Enter  data  using  DEI  Sequence  (Alternate  entry i  store  n  in  R, 

2  2 
sx  in  R^,  m  in  and  s^  in  RQ^) 

2.  Store  either  ( l+'Y)/2  in  Rq^  OR  F^+^^2(n-l  ,m-l)  in  R^ 

AND 

F(l+7')/2^m”1,ri"1^  ln  R13 

3*  Press  1  is  displayed,  then 

press  [x^t |  u  is  displayed 
C.I.  For  EXPONENTIAL  \X/\Y  =  flY/flx 

1.  Enter  data  using  DEI  Sequence  (Alternate  entry:  store  n  in  R1 
x  in  R^g,  m  in  RQ^  and  y  in  RQg) 

2.  Store  either  ( l+'y)/2  in  RQ^  OR  Fj-1+-yy2(2m,2n)  in  R^ 

AND 

F(l+7')/2^2n,2m)  ^  R13 

3.  Press  [  s'  I  1  is  displayed,  then 
press  fx£t|  u  is  displayed 

PROGRAM  2  ADDITIONAL  CAPABILITIES: 


Inverse  Normal  CDF  Approximation 

1.  Store  p  in  R^  (p  =»  .5) 

2 .  Press  □D  to  display  z^ 


Inverse  Student's  t  CDF  Approximation 


1. 

Store  p  in  RQ9  (p 

>  .5)  and  v  in  R^  (v  ^ 

1) 

2. 

Press  SBR  ; 

RCL 

to  display  t  (v) 

Inverse 

F  CDF  Approximation 

1. 

Store  p  in  RQ9  (p 

>  .5),  vt  in  R2q  and  v£ 

in  R2i 

2. 

Press  SBR  STO 

to  display  Fp(vi,V2) 

(V^V1) 
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PROGRAM 

2  LABELS  USED 

A 

A’ 

GTO 

EE 

B 

3' 

RST 

G 

= 

SBR 

D 

D' 

STO 

r-i 

2m 

E* 

RCL 

PROGRAM 

2  STORAGE  REGISTER 

CONTENTS : 

00 

clear 

15 

n 

01 

nx  or 

16 

used 

02 

my  or  £3^  2 

17 

s/  or 

Ox2/n 

03 

non 

18 

X 

04 

m  or  JXi 

19 

v  or  (n 

-1) 

05 

IX2 

20 

V1 

06 

Tvi 

21 

V2 

07 

sy2  or  (JY2/m 

22 

clear 

08 

y 

23 

used 

09 

p  or  (l+7)/2 

24 

used 

10 

x-y 

25 

clear 

11 

GDF  value 

26 

clear 

12 

8aV“" 

27 

used 

13 

CDF  value 

28 

used 

14 

used 

29 

used 

PROGRAM  2  SAMPLE  PROBLEMS  ; 

1.  In  a  survey  of  400  people  from  one  city  188  preferred  Brand  A  soap 
to  all  others;  and  in  a  sample  of  500  people  from  another  city  210 
preferred  the  same  product .  Find  a  95%  confidence  interval  for 

PX_PY  ’ 

SOLUTION; 

(1)  Press  then  store  nx  =  188  in  Rq  ,  my  =  210  in 

n  =  400  in  RQ^  and  m  =  500  in  RQ^ 

(2)  Store  (l+y)/2  =  .975  in  RQ9  OR  z_9?5  =  1 .960  in  R^ 

(3)  Press  m  then  |x£t |  to  display 

1  =  -.0153123448  1  =  -.0152991877 

OR 

u  =  .1153123448  u  =  .1152991877 

2.  In  a  survey  of  dieting  effects,  ten  women  were  selected,  weighed 
and  placed  on  a  diet  for  two  weeks.  At  the  end  of  that  time  they 
were  reweighed.  The  results  are  listed  below.  Find  a  9 9%  confidence 
interval  for  the  difference  In  means .  Also  find  a  95%  upper 
confidence  bound  for  the  difference  in  means . 

Before  119  122  136  130  129  136  134  133  119  115 

After  114  119  1>  126  119  137  124  127  119  107 

SOLUTION  Part  a. 

(1)  Enter  data  using  DEP  Sequence,  9  is  displayed 

(2)  Store  (1+70/2  =  .995  in  RQ9  OR  t  .995(9)  =  3-250  in  Ru 

(3)  Press  m  then  [xgt |  to  display 

1  =  .7054594846  1  =  .7328694249 

OR 

u  =  8.694540515  u  =  8.667130575 
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SOLUTION  Part  b. 


(1)  Enter  data  using  DE?  Sequence,  9  is  displayed 

(2)  Store  J=  .95  in  RQ9  OR  t  <95(9)  =  1-833  in  Ru 

(3)  Press  m  then  j  x<t  ]  to  display 

ignore  1  ignore  1 

OR 

u  =  6.92957002  u  =  6. 937461 644 

3-  Suppose  x  =  10,  y  =  5.  n  =  25,  and  s^  =  10;  find  a  90%  C.I.  for 

Mx'My  ‘ 

SOLUTION! 

(1)  Press  m  then  s  ore  x-y  =  5  in  R^q,  sd/Vr  =  2  in 
and  (n-1)  =  24  in  R^ 

(2)  Store  (l+7)/2  in  R^  OR  1^(24)  =  1.711  in  R^ 

(3)  Press  m  then  fxVt|  to  display 

1  =  1.581803572  1  =  1.578 

OR 

u  =  8.418196428  u  =  8.422 

4.  Two  small  classes  used  different,  but  equivalent  methods  on  a 

common  exam .  Their  scores  are  listed  below .  Find  a  90%  confidence 
interval  and  a  90%  lower  bound  for  the  difference  in  means. 

X:  82  87  91  54  97  76  64  98  92  57  80  53  64 

Y:  91  62  94  92  87  79  86  75  90  73  83  93  65  89  68  52 

SOLUTION;  Part  a. 

(1)  Enter  data  using  DEI  Sequence,  v  =  27  is  displayed 

(2)  Store  (l+y)/2  =  .95  in  RQ9  OR  t>95(27)  =  1-703  in  Rn 

(3)  Press  m  then  [x%t |  to  display 

1  =  -12.62699316  1  =  -12.63357801 

OR 

u  =  5.828916241  u  =  5.835501092 
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SOLUTION  Part  b. 

(1)  Enter  data  using  DEI  Sequence,  v  =  27  is  displayed 

(2)  Store  J  =  .90  in  OR  t  (27)  =  1-314  in  Rn 

(3)  Press  m  to  display 

1  =  -10.49428744  1  =  -10.52422048 

OR 

ignore  u  ignore  u 

5.  Suppose  n  =  16,  m  =  15,  x  =  14.3,  y  =  12.5.  ^(x^-x)2  =  67*2  and 
J](yi-y)2  =  103.7  *  Find  a  99%  upper  confidence  bound  for  the 
difference  in  means . 


SOLUTIONS 

(1)  Press  m  then  store  n  =  16  in  R^>  m  =  15  in  Rq^,  x- 
1.8  in  R10,  ^(xi-x)2  =  67.2  in  Rl6,  J^-y)2  =  103-7 
in  R26  and  (n+ra-2)  =  29  in 

(2)  Store  J  =  .99  in  RQ9  OR  t  (29)  =  2.462  in  Rn 

(3)  Press  H3  then  jx^t I  to  display 

ignore  1  ignore  1 

OR 

u  =  3-959764742  u  =  3-948005095 

6.  Suppose  n  =  32,  m  =  10,  x  =  13,  y  =  15,  -  16  and  (J^2  =  25; 

find  a  95%  confidence  interval  for  Mx'My* 

SOLUTION s 

(1)  Press  lD  then  store  x-y  =  2  in  and 
/(CTx2/n)  +  (CTY2 A)  =  ifT in  r12 

(2)  store  (l+y)/2  =  .975  in  RQ9  OR  z<9?3  =  I.96O  in  Ru 

(3)  Press  m  then  |x%t|  to  display 

1  =  -1.395303599  1  =  -1-394819583 

OR 

u  =  5-395503599  u  =  5-394819583 


65 


M 


^  I 


1 


7-  Given  the  following  sets  of  data  with  known  variance,  estimate  the 
difference  in  means  with  a  90%  confidence  interval . 

<JX2  =  3,7  and  g"Y2  =2.5 
X:  3.6  4.7  2.4  1.7  6.2  7-2  3.6 

Y:  5-4  7.6  4.3  6.5  8.7  4.5  6.1  8.6 

SOLUTION! 

(1)  Enter  data  using  DEI  Sequence  and  store  0^/n  =  3-7/7  in 
R1?  and  (JY2/m  =  2.5/8  in  RQ? 

(2)  Store  (1+7 *)/2  =  .93  in  RQ9  OR  z  ^  =  1 .643  in  R^ 

(3)  Press  m  then  | xVt|  to  display 

1  =  -3-771322331  1  =  -3.77112862 

OR 

u  =  -.7536774686  u  =  -.7538713802 

8 .  Rind  a  99^  confidence  interval  for  £T2/  CTV 2  using  the  following 
data . 

X:  82  87  91  54  97  76  64  98  92  57  80  53  64 

Y:  91  62  94  92  87  79  86  75  90  73  83  93  65  89  68  32 

SOLUTION : 

(1)  Enter  data  using  DEI  Sequence 

(2)  Store  (i+y)/2  =  .995  in  R^  OR  F.995(12,15)=4.25  in  Rn 

AND 

F_995(15,12)  =  4.72  in  R13 

(3)  Press  m  then  |~x%t|  to  display 

1  =  .3864689829  1  =  .3875258IO8 

OR 

u  =  7-818933411  u  =  7-773767765 
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.  Find  a  confidence  interval  for  (Jx  /  Qy  given  n  =  3,  m  =  11, 
s  2  =  4  and  s  2  =3-6  . 

A  1 

SOLUTION: 

(1)  Press  1  D  j  th®n  store  n  =  8  in  R^,  m  =  11  in 
s„2  =  4  in  R  and  s ^  =  3*6  in  RQ7 

(2)  Store  (l+7)/2  =  -975  in  RQ9  OR  F<975(7,10)  =  3-95  in  R^ 

AND 

Fi9?5(10,7)  =  4.76  in  R^ 

(3)  Press  mu  then  |x%t)  to  display 

1  =  .280241493  1  =  .2812939522 

OR 

u  =  5-35837828  u  =  5.288888889 

0 .  Find  a  9 0 %  confidence  interval  for  M/M  ^  given  the  following 
data: 

X:  17  6  12  14  3  12  15 

Y:  7  15  4  6  21  13 

SOLUTION: 

(1)  Enter  data  using  DEI  Sequence 

(2)  Store  (l+7)/2  =  .95  in  RQ9  OR  F>95(12,14)  =  2.55  in  Ru 

AND 

F>95(14,12)  =  2.64  in  R^ 

(3)  Press  m  then  [x%t |  to  display 

1  =  .3688954123  1  =  .3822288409 

OR 

u  =  2.472430016  u  =  2.5731o4557 


11.  Two  types  of  electric  bulbs  axe  observed  as  to  length  of  life, 
yeilding  the  following  results.  Are  the  means  significantly 
different?  (C(  -  -1) 

Type  1  n  =  46  x  =  1070 
Type  2  m  =  64  y  =  1041 

SOLUTIONS 

(1)  Press  m  then  store  n  =  46  in  R^,  m  =  64  in  RQ^, 
x  =  1070  in  R^g  and  y  =  1041  in  RQg 

(2)  Store  (l+y)/2  =  .95  in  RQ9  OR  F  .^(98,123)  =  1.33  in  Rn 

AND 

Fi95(128,98)  =  1.40  in  R 

(3)  Press  □D  then  |x£t]  to  display 

1  =  .7102715868  1  =  .7315016513 

OR 

u  =  1.345613548  u  =  1.362056075 

Thus  with  90%  confidence  rhe  means  axe  not  different . 


PROGRAM  3  USER  GUISE  -  One -Population  Hypothesis  Tests 
INTRODUCTION;  The  purpose  of  this  program  is  to  perform  one-tailed  or 
two-tailed  hypothesis  tests  for  the  following  one-population  situations; 
NORMAL  with  (J2  known 

NORMAL  fj LQ  with  (J2  unknown 

BERNOULLI  p 

o 

NORMAL  (J2  with  f±  known 

NORMAL  <J2  with  fl  unknown 

EXPONENTIAL  (jLo  =  l/^Q 
POISSON  \Q 

GENERAL  PROCEDURES; 

1 .  This  program  requires  the  Applied  Statistics  Module .  After 
reading  all  three  card  edges,  press  EH  to  repartition  ( 71 9 .29). 

2.  For  each  new  problem  press  CE1  to  clear  all  registers  and  to 
prepare  for  the  following  parameter  entries.  Enter  Qq  and  press  R/s  , 
enter  CL  and  press  r/s  ,  then  enter  either  0,  -1  or  +1  for  the  desired 
alternate  hypothesis,  as  in  the  table  below,  and  press  r/s  • 


0 

for 

V 

e  /  eQ, 

Two-tailed  test 

-1 

for 

H1‘ 

e  -  e0- 

Lower-tailed  test 

+1 

for 

V 

e  -  ©o* 

Upper-tailed  test 

3.  For  data  entry  press  m  followed  by  data  point  x^,  [r/s! , 

xi+l’  ft/S  >  e"ta  for  each  x^^  (i=l  ,2 . n)  until  all  data  points  have 

been  entered.  Mistakes  in  data  entry  should  be  corrected  immediately  by 
reentering  the  unwanted  data  point  and  pressing  INV  [2nd|  E3  ,  then 
enter  the  correct  data  point  and  press  R/s  f r/ s| .  Alternate  data  entry 
using  summary  statistics  is  detailed  in  each  subroutine . 
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4.  At  the  conclusion  of  each  test  a  1  for  Reject  or  a  0  for 
Accept  is  displayed.  The  significance  level  for  each  test  is  usually  in 
the  T-register  (R  ) ;  however,  for  two-tailed  unsymmetric  tests  the  test 
statistic  stored  in  R^  is  used  as  follows :  if  a  0  is  displayed  then 
the  test  must  have  failed  to  reject  in  both  tails.  If  a  1  is  displayed 
then  the  test  was  upper-tailed  if  R^  is  greater  than  the  median  and 
lower-tailed  if  R^  is  less  than  the  median. 

PROGRAM  3  SPECIFIC  PROCEDURES : 

Tests  for  NORMAL  f_lQ  with  (J^  known 

1 .  Press  mn  ,  enter  jJLQ,  press  |r/ S j  ,  enter  a,  press  }r/s| , 
enter  0,  -1  or  *1  and  press  R/s  . 

2.  Enter  data  using  and  store  (J  in  RQ^  (Alternate  entry: 

store  n  in  R^,  x  in  R^q  and  (J  in  R^)  . 

3-  Press  [  A  |  either  1  (reject)  or  0  (accept)  is  displayed, 


press 


to  display  significance  level. 


Tests  for  NORMAL  fJ,Q  with  (J^  unknown 

1 .  Press  H’  enter  press  [r/s) »  enter  (X,  press  |r/s| , 
enter  0,  -1  or  +1  and  press  R/s  . 

2.  Enter  data  using  EJ  (Alternate  entry:  store  n  in  Rq^»  x  in 
and  Sj(  in  R^)  . 

3.  CD  either  1  (reject)  or  0  (accept)  is  displayed,  then 


pre  ss 


to  display  significance  level. 


Tests  For  Bernoulli  p 

- - 

1 .  Press  enter  p  ,  press  (r/s | ,  enter  0( ,  press  |r/s|  , 

enter  0,  -1  or  +1  and  press  R/s  . 

2.  Store  nx  in  R^  and  n  in  Rq^. 

3 .  Press  m  either  1  (reject)  or  0  (accept)  is  displayed,  then 
press  |xlt j  to  display  significance  level. 

Tests  For  NORMAL  (J2  with  fj.  known 

1 .  Press  cm,  enter  (J2 ,  press  (r/s]  ,  enter  d  ,  press  ffl/sj  , 
enter  0,  -1  or  +1  and  press  |r/s|  • 

2.  Enter  data  using  cm  and  store  fJi  in  R^q  (Alternate  entry: 
store  n  in  R ,  JV  in  RQ1 ,  Jx^2  in  R02  and  j. I  in  R^Q)  . 

3*  Press  f  D ' |  either  1  (reject)  or  0  (accept)  is  displayed,  then 
press  |x£t |  to  display  significance  level . 

Tests  Far  NORMAL  (J2  with  jj.  unknown 

1 .  Press  [  E^] ,  enter  Q mQ2  ,  press  r/s  ,  enter  (X ,  press  P./s  i 
enter  0,  -1  or  +1  and  press  R/s  • 

2.  Enter  data  using  m  (Alternate  entry:  store  n  in  RQ^  and 

(n-l)sX  R10^  * 

3 ■  Press  nn  either  1  (reject)  or  0  (accept)  is  displayed,  then 
press  [x%t|  to  display  significance  level. 

Tests  For  EXPONENTIAL  j±Q  =  l/^Q 

1 .  Press  cm.  enter  jj.  0>  press  [r/s]  »  enter  Qi,  press  | R/s |  , 
enter  0,  -1  or  +1  and  press  R/s  . 

2.  Enter  data  using  )  D  ]( Alternate:  store  n  in  R^^  and  *  in  R10)  • 

3.  Press  [  E  ]  either  1  (reject)  or  0  (accept)  is  displayed,  then 
press  \x$X  ]  to  display  significance  level . 
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Tests  For  POISSON  \ 

1 .  Press  [ID  ,  enter  \  ,  press  ]r/s] ,  enter  Oi  ,  press  |r/s|  , 
enter  0,  -1  or  +1  and  press  R/s  . 

2.  Enter  data  using  [~*D~]( Alternate  entry j  store  n  in  Rq^ 
nx  in  Rq1) . 

3 .  Press  cm  either  1  (reject)  or  0  (accept)  is  displayed,  then 
press  |x<t |  to  display  significance  level. 

PROGRAM  3  LABELS  USED 


A 

A' 

TAN 

LNX 

3 

3' 

SIN 

LOG 

G 

C' 

COS 

INV 

D 

D' 

CLR 

FIX 

E 

E' 

GRD 

INT 

X 

+ 

SUM 

ENG 

_•  _ 

• 

- 

GE 

PROGRAM  3  STORAGE  REGISTER  CONTENTS 

CO  0,-1  or  +1  10  x  or  fj.  or  (n-l)s^ 

01  nx  or  £x^  11  a 

02  Yx.Z  12  test  statistic 


02 

03  n 
04  clear 
05  clear 

06  degrees  of  freedom 

0?  (J 

08  flQ  or  CT02  or  PQ  or  \ , 
09  clear 


13  clear 

14  sx 

15-29  used 


PROGRAM  3  SAMPLE  PROBLEMS  : 

1 .  According  to  an  early  encyclopedia  the  average  rainfall  in  a  city  is 
30.1  inches.  Rainfall  during  the  past  five  years  has  been: 

30.5  27.4  35-1  32.6  25.9 


Assuming  a  standard  deviation  of  .2  inches,  has  the  average  changed? 

SOLUTION:  Hq:  f±=  30.1  ,  Hj :  30.1  ,  CL  =  -05 

(1)  Press  ID>  30.1,  |r/s! ,  .05,  [r/s|,  0,  [r/s! . 


•  (2)  Enter  data  using  [~D  |,  store  (J  =  .2  in  RQ7 

(3)  Press  CD  1  is  displayed  (reject  H  ), 
fx£tj  si  =  .0253472347  . 

2.  A  new  fad  diet  was  tried  out  on  15  subjects  and  the  weight  losses 
after  one  week  were: 


2.07  2.34  1.97  1.85  1.84  2.23  2.15  1.89 

1.93  1-99  1.86  1.90  2.09  2.16  2.04 

An  advertisement  claims  that  the  weight  loss  after  one  week  on  the 
diet  is  at  least  two  pounds.  Do  the  data  support  the  claim ?(0(=  .05) 
SOLUTION:  Hq:  ^  22  ,  ^ :  fl  -e  2 

(1)  Press  mu  ,  2,  [r7s],  .05,  [r/s],  1, 

(2)  Enter  data  using  0. 

(3)  Press  DO  0  is  displayed  (accept  H  ) , 

\xtt\  si  =  .302539836  . 

3.  It  is  claimed  that  a  certain  drug  will  lower  temperature  within  ten 
minutes.  Five  subjects  having  normal  temperatures  of  98.6  were  given 
the  drug  and  ten  minutes  later  their  temperatures  were  recorded  in 
summary  form  as  follows:  ]Tx  =  491  ,  ]Tx2  =  48219  •  Test  the  claim 
at  the  Ol  =  .05  level. 


r/s 
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SOLUTION;  Hqj  ^2  98.6  ,  H^/is-98.6 

(1)  Press  GD  ,  98.6,  [r/s],  .05,  |r/s[,  -1,  |R/sl ■ 

(2)  Store  n  =  5  in  R^.  x  =  491/5  in  R10  and  in  R^ 

(3)  Press  CD  0  is  displayed  (accept  H0) , 

|xVt]  si  =  .1352267548  . 

4.  A  student  claims  he  can  always  answer  more  than  half  of  the  items  on 
a  true-false  exam  correctly,  regardless  of  the  topic.  You  devise  a 
20  question  exam  on  a  subject  of  which  he  knows  nothing  and  he 
answers  12  of  the  test  items  correctly.  Would  you  conclude  he  has 
extraordinary  powers?  (  0(  =  .05) 

SOLUTION;  H q;  /!=  .5  ,  Hj  s  [l  =*  -5 

(1)  Press  H  ■  .5,  [575],  *05,  [575],  i.  [r/s]  • 

(2)  Store  nx  =  12  in  RQ^  and  n  =  20  in  R„. 

(3)  Press  m  0  is  displayed  (accept  Hq) , 

[xS.t  |  si  =  .2517223358  . 

5.  Of  694  respondents,  3^9  were  in  favor  of  more  dissemination  of  birth 
control  Information.  Is  it  safe  to  conclude  that  more  than  half  of 
the  population  agrees  with  this  position?  (Oi=  .05) 

SOLUTION;  Hq:  fj.  S  .5  ,  H 1*  JJL  -5 

(1)  Press  GO  .  .5,  [r/s],  .05,  [r/s],  1,  US]  • 

(2)  Store  nx  =  369  in  Rqi  a^d  n  =  694  in  Rq^ 

(3)  Press  CD  1  is  displayed  (reject  HQ) , 

|xStl  si  =  .0474381802  . 


6.  A  soup  can  filling  machine  is  supposed  to  fill  each  can  with  ten 
ounces  of  clear  broth,  with  a  variance  of  .01 .  A  change  in 
variability  in  either  direction  is  undesirable .  A  random  sample  of 
20  cans  yeilded  £x2  =  2107.?  and  Jx  =  205-3  •  Is  the  machine 
working  within  limits?  (OC=  *05) 

SOLUTION  t  Hq:  <J 2  =  .01  ,  ^  s  <J2  t  .01 

(1)  Press  E1 \  ,  .01,  |r/s[,  .05,  |r/s| ,  0,  (r/s[  . 

(2)  Store  n  =  20  in  R^,  Jx  =  205-3  in  RQ1 ,  £x2  =  210?.?  in 

R02  and  M  =  10  in  R10 

(3)  Press  [  D ' |  1  is  displayed  (reject  Hq) . 

7.  A  nail  machine  is  supposed  to  manfacture  1-inch  nails  with  a 
standard  deviation  of  .025  inches.  A  random  sample  of  3°  nails 
yeilded  a  sample  value  for  sY  of  .03  inches.  Does  this  apparent 

A 

increase  warrant  shutting  the  machine  down?  (  =  *°5) 

SOLUTION;  HQ:  (J2  58  (.025)2  ,  Hj.  ;  (J  2  s-  (  .025)2 

(1)  Press  tun .  .000625,  [r/s[ ,  .05,  [r/s|,  1,  |r/s! . 

(2)  Store  n  =  30  in  Rq^  and  (n-l)s2  =  29( .O3)2  in  R1Q 

(3)  Press  El  0  is  displayed  (accept  Hq) 

fcjt]  si  =  .059 

8 .  A  certain  type  of  expensive  electrical  gear  is  supposed  to  have  a 
mean  life  of  1000  hours .  The  manfacturer  is  concerned  if  the  mean 
departs  in  either  direction  from  1000.  Five  components  were  tested 
and  they  had  the  following  burnout  times; 

1075  1085  1060  998  995 

Is  the  mean  still  1000?  ( CL  =  -05) 


SOLUTION!  Hq!  11=  1000  ,  ^ :  fl  4  1000 

(1)  Press  GD  ,  1000,  [r/s! ,  .05,  [r/s| ,  o,  [r/s]  ■ 

(2)  Enter  data  using  tH 

(3)  Press  H  0  is  displayed  (accept  H  ) 

9  •  Over  a  period  of  years  there  had  been  an  average  of  14  accidents  per 
year  in  a  certain  city.  This  year  the  monthly  totals  were  as 
follows!  10221  1301012. 

Does  this  data  agree  with  the  theory  that  the  number  of  accidents 
per  month  follows  a  poisson  distribution  with  \=  1.166??  (Oi~  .05) 

SOLUTION!  Hq  *•  \=  1.1667  ,  VX*  1.1667 

(1)  Press  nn  ,  1.1667,  [rTs],  .05,  [r/s ) ,  0,  [3/s I  . 

(2)  Enter  data  using  m 

(3)  Press 


GO  0  is  displayed  (accept  H  ) 
|^jt]  si  =  .5703943172 


PROGRAM  4  USER  GUIDE  -  Two-Population  Hypothesis  Tests 
INTRODUCTION i  The  purpose  of  this  program  is  to  perform  one-tailed  and 
two-tailed  hypothesis  tests  for  the  following  two-population  situations. 
BERNOULLI  %=PY  for  larSe  sample  sizes 

NORMAL  Mx=/^Y  f°r  X  and  Y  Paired 

NORMAL  Px=/Xy  ~or  X  and  *  independent  with  (Jy^= 

NORMAL  /^X=M Y  for  X  and  Y  independent  with  Q J™0*11 

NORMAL  Mx=My  Yor  *  and  independent  with  (Jx2  /  (J^2 

NORMAL  CTv2=C7V2  ^or  x  ar-d  Y  independent 

NORMAL  p  =  0  for  X  and  Y  paired 

EXPONENTIAL  X^Xy  Yar  Y  and  ^  independent 

GENERAL  PROCEDURES; 

1.  This  program  requires  the  Applied  Statistics  Module.  After 
reading  all  three  card  edges,  press  E  to  repartition  (719.29). 

2.  For  each  new  problem  press  E  to  clear  all  registers  and  to 
prepare  for  the  following  parameter  entries.  Enter  Qi  and  press  R/s 
then  enter  either  0,  -1  or  +1  for  the  desired  alternate  hypothesis,  as 
in  the  table  below,  and  press  r/s  . 

0  for  :  Q ^  ^  ,  Two-tailed  test 

-1  for  s  «e  ,  Lower-tailed  test 

+1  for  i  s»  Qy  ,  Upper-tailed  test 

3-  For  independent  data  entry  (DEI  Sequence)  press  S  followed 
by  data  point  x^,  r/s  ,  x^  ,  (r/ s|  ,  etc.  for  each  x^  (i=l  ,2,  .  .  .  ,n)  . 
When  all  x^'s  have  been  entered  press  J~~D~] again,  followed  by  data  point 
yi’  R/s  *  yi+l’  R/s  *  etc.  for  each  (i=l,2, . . .,m) .  when  all  of  the 
yj_’s  have  been  entered  press  |  D  j.  Mistakes  in  data  entry  should  be 
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corrected,  immediately  by  reentering  the  unwanted  data  point  and  pressing 
INV 1  [2nd]  S3  ,  then  enter  the  correct  data  point  and  press  [R/S| 
r/s  .  Alternate  data  entry  is  detailed  in  each  subroutine. 

4.  For  paired  data  entry  (DEP  Sequence)  press  m  followed  by  x^, 
|xit  j ,  y.  |R/s| ,  x,  ,  „  ,  fxit|  ,  y.M,  |r/s|  ,  etc .  when  all  of  the  paired 


.  xi+1 ,  [x|tj ,  yi+1 , 


data  points  have  been  entered  press  [ SER  j  [ RST |  .  Mistakes  and  alternate 
entry  are  as  above . 

5.  Each  routine  displays  either  a  1  for  Reject  or  a  0  for 
Accept .  The  significance  level  is  always  in  the  T-register  (RT) . 

PROGRAM  4  SPECIFIC  PROCEDURES! 

Tests  for  BERNOULLI  PX"PY  for  larSe  sample  sizes 

1 .  Press  HD  ,  enter  a,  press  |r/s[ ,  enter  C ,  -1  or  +1  ,  and 
press  |r/s|  . 

2.  Store  n  in  R^,  nx  in  R^^,  m  in  RQ^  and  my  in  RQ1  . 

3 .  Press  m  either  1  (reject)  or  0  (accept)  is  displayed,  then 
press  EH  to  display  significance  level . 

Tests  for  NORMAL-  for  X  and  Y  paired 

1 .  Press  m  ,  enter  a,  press  [575],  enter  0,  -1  or  +1  and  then 
press  |r/ S I  • 

2.  Enter  data  using  DE?  Sequence  (Alternate  entry:  store  (n-1)  in 
R25  and  in  R10^ 

3*  Press  □3  either  1  (reject)  or  0  (accept)  is  displayed,  then 


press  x£t  to  display  significance  level. 


Tests  for  NORMAL  =  j±Y  with  (j*  =  <jY2  =  (J 2 

1 .  Press  El  ,  enter  oc,  press  |r/s[ ,  enter  0  ,  -1  or  +1  and  then 
press  |'r/S  | . 

2.  Enter  data  using  DEI  Sequence  (Alternate  entry  j  store  n  in  R^, 
n  in  RQy  x  in  R^^,  y  in  Rq8,  sx2  in  R^  and  sy2  in  R^) 

3  •  Press  E  either  1  (reject)  or  0  (accept)  is  displayed,  then 
press  S3  to  display  significance  level . 

Tests  for  NORMAL  fj for  q-^2  and  known 

1 .  Press  E  ,  enter  (X,  press  |R/s| ,  enter  0,  -1  or  +1  and  then 
press  [r/s{  . 

2.  Enter  data  using  DEI  Sequence  and  store  (J^  in  R^  and  (J^2  in 
Rq^  (Alternate  entryi  store  n  in  R^,  m  in  RQ^,  x  in  R^,  y  in 

^08’  C?"x  ln  ^09  and  ^Y2  in 

3 .  Press  m  either  1  (reject)  or  0  (accept)  is  displayed,  then 
press  | x^t |  to  display  significance  level. 

Tests  for  NORMAL  =  ££y  for  O /  (Jy2 

1 .  Press  E  ,  enter  a>  press  [r/s| ,  enter  0,  -1  or  +1  and  then 
press  f R/S [ . 

2.  Enter  data  using  DEI  Sequence  (Alternate  entry:  store  n  in 
m  in  Rq^ i  x  in  Rli+,  y  in  RQg,  Sy2  in  R^  and  Sy2  in  R^) . 

3 .  Press  E  ,  1  (reject)  or  0  (accept)  is  displayed,  then 
press  | x^t|  to  display  significance  level . 


Tests  for  NORMAL  (JZ  =  (J.Z 


1 .  Press  □3  ,  enter  a,  press  | r/s { ,  enter  0,  -1  or  +1  and  then 


press 


2.  Enter  data  using  DEI  Sequence  (Alternate  entry  j  store  n  in  R^, 

2  2 

di  in  Rq^i  in  Rq^  2-nd.  Sy  in  Rq^)  • 

3 .  Press  nn  either  1  (reject)  or  0  (accept)  is  displayed,  then 
press  |x£t j  to  display  significance  level. 

Tests  for  NORMAL  f)=  0 

1 .  Press  GD  ,  enter  CL,  press  [r/ S  j ,  enter  0 ,  -1  or  +1  and  then 


press  | R/S |  . 

2.  Enter  data  using  DEP  Sequence  (no  alternate  entry). 

3 .  Press  E  either  1  (reject)  or  0  (accept)  is  displayed,  then 
press  | x%t [  to  display  significance  level. 

Tests  for  EXPONENTIAL  Xx  =  Ay 

1 .  Press  GD  ,  enter  a,  press  [r/s | ,  enter  0,  -1  or  +1  and  then 
press  [r/s[ . 

2.  Enter  data  using  DEI  Sequence  (Alternate  entry;  store  n  in  R^, 
m  in  RQy  x  in  R^  and  y  ir.  R^g)  . 

3 .  Press  □Q  either  1  (reject)  or  0  (accept)  is  displayed,  then 
press  [xStj  to  display  significance  level. 

PROGRAM  4  LABELS  USED 


A 

A’ 

SIN 

RST 

GTO 

B 

B' 

COS 

CLR 

X^ 

G 

G' 

DEG 

RCL 

+ 

D 

D’ 

INV 

STO 

- 

E 

E' 

GRD 

EXC 

PROGRAM  4  STORAGE  REGISTER 

CONTENTS; 

00 

0,  -1  or  +1 

07 

2  ? 
sY  or  cr/ 

14  x 

01 

my  or  ]Ty 

08 

y 

15  n 

02 

I*2 

09 

s/  or  cr/ 

16-24  used 

03 

m  or  n 

10 

25  (n-1) 

04 

nx  or  £x 

11 

a 

26  -  29  used 

05 

I*2 

12 

used 

06 

13 

used 

PROGRAM 

4  SAMPLE  PROBLEMS ; 

1.  According  to  last  year.s  statistics,  70%  of  the  population  were  in 
favor  of  stricter  smog  control  laws .  This  year  only  of  1000 
respondents  favored  stricter  laws .  Does  this  represent  a  significant 
decrease?  (  Q£=  -oi) 

SOLUTION;  Hq;  Px  =  PY  •  % .  Px  Py 

(1)  Press  E  ,  .oi,  [i/s] ,  +i, 

(2)  Store  n  =  1000  in  R^f  nx  =  700  in  Rq^,  m  =  1000  in  RQ^ 
and  my  =  6 50  in  RQ^ . 

(3)  Press  m  1  is  displayed  (reject  H  ) 

E3  si  =  .0084920904 

2.  Ten  plots  are  split  and  half  planted  with  variety  a  ar.d  half  with 
variety  B .  The  yields  are  shown  below ,  is  there  a  difference  in  mean 
yield  between  the  two  varieties?  (a=  -io) 


A;  49 

58 

53 

60 

45 

49 

66 

55 

44 

52 

B;  47 

57 

49 

57 

44 

44 

67 

52 

42 

53 

sou™.  V  Ma  “Mb  -  K1  Ma  I'  Mb 
(1)  Press  E  ,  .io,  |r/s]  ,  o,  [57s] . 


R/S 
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(2)  Enter  data  using  DEP  Sequence. 

(3)  Press  cm  1  is  displayed  (reject  H  ) 

1 x^t |  si  =  .0137674905 

3 •  A  cigarette  manufacturer  tests  tobacco  from  two  different  brands  for 
nicotine  content  and  obtains  the  following  (in  milligrams)  . 

Ai  24  26  25  22  23 

B:  27  28  25  29  26 

Do  these  results  indicate  there  is  a  difference  in  mean  nicotine 
content  for  the  two  brands?  (a=  -05) 

SOLUHON:  Hq:  =  MS  ’  H1  ‘  M A  ^3 

(1)  Press  s  -  .05,  fi/s|,  0,  fj/s] . 

(2)  Enter  data  using  DEI  Sequence. 

(3)  Press  m  1  is  displayed  (reject  H  ) 

fx^t]  si  =  .0170716812 

4.  Two  analysts  make  fifty  independent  determinations  of  the  melting 
point  of  a  certain  chemical .  The  sample  mean  and  variance  of  the 
data  found  by  analyst  I  are  respectively,  73-6  and  10  while  the 
sample  mean  and  variance  found  by  analyst  II  are,  respectively, 

72 .4  and  8  .  It  is  argued  that  there  is  a  tendency  for  analyst  I 
to  get  higher  results.  What  is  your  conclusion?  {0(=  .05) 

SOLUTIOT.  V  /4'My’ 

(1)  Press  cm  .  .05,  [r7s],  1,  [hTH  • 

(2)  Store  n  =  50  in  R^,  m  =  50  in  R^,  x  =  73-6  in  R^, 

Y  -  72.4  in  Ro8,  sx2  =  10  in  R09  and  sy2  =  8  in  R0?. 

(3)  Press 
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m  1  is  displayed  (reject  H  ) 
|x%t]  si  =  .022750062 


s 


5 


r 


*  «. 


1 


We  are  interested  in  replacing  wire  B  with  wire  A  if  the  resistance 
per  unit  length  is  not  significantly  decreased.  The  results  of 
twenty  tests  on  each  wire  are  presented  below.  If  we  know  that  the 
standard  deviation  of  the  two  testing  procedures  are  both  .001?  ohms. 
What  is  your  recommendation?  (  0^  =  .01) 

A 


.051 

.051 

.049 

.048 

.049 

.049  .049  .053 

.053 

.049 

.049 

.047 

.048 

.049 

.049 

.051  .051  .050 

.047 

.050 

.054 

.052 

.051 

.049 

.051 

.057  .054  .052 

.050 

.052 

.051 

.051 

•  055 

.049 

.052 

.051  .051  .052 

.052 

.048 

SOLUTION  1 

V 

Ma  =  (1b  • 

H1  ‘  Ma  >  * 

(1) 

Press 

L*i 

,  .01, 

R/sI 

|,  1,  |  r/s  . 

(2) 

p 

Enter  data  using  DEI  Sequence  and  store  (J^~ 

=  ( .0017)' 

ln  R09 

:  and 

(Jy2  = 

(  .0017)2  in  Rq? 

(3) 

Press 

d 

0  Is  displayed  (accept  Hq) 

1  si  =  .500837784 

6.  Twenty  plots  of  ground  were  planted  with  corn.  Ten  plots  (Y) 

contained  a  special  treatment .  The  variances  were  tested  and  found 
to  be  unequal .  Is  there  a  significant  difference  between  the  yields 
at  the  5%  level? 

x  =  6.1  y  =  5-78 


sx2  =  .13556 


sy2  =  .02844 


1 

SOLUTION: 

rH 

K 

II 

:£ 

O 

?! 

(1) 

Press 

E’  ,  .05,  R/S  ,  0, 

- 1 

4 

(2) 

Store  n  =  10  in  R^,  m  =  10 

< 

i . 

•1 

(3) 

in  R0g, 

Press 

,  s/  =  .13556  in  R09 
{  G '  |  1  is  displayed 

03’ 

V 


"14  ’ 
R07 ' 


[  xttl  si  =  .026650012 
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7.  Suppose  that  two  samples  of  ten  and  sixteen  observations  have, 
respectively,  variances  of  .3888  and  2.25  •  At  the  5%  significance 

level  would  you  accept  the  hypothesis  that  (J ^2  <  (J£2  ? 

SOLUTION!  Ho;  (Jx2s  crY2  *  V  0*x2=*C7y2  • 

(1)  Press  EZ1  •  .05,  |r/s  1 ,  -1,  |r/s|  . 

(2)  Store  n  =  10  in  R  ,  m  =  16  in  Rq  ,  sx2  =  .3888  in  R  and 
sy2  =  2.25  in  RQ?. 

(3)  [  D' |  0  is  displayed  (accept  H  ) 

[x*t]  si  =  .99^190704  . 

8.  The  following  data  measurements  of  students’  ability  in  an  IQ  test 
are  paired  with  scores  in  an  acheivement  test.  Test  for  p  =  0. 

Xi  105  95  125  92  120  10?  121  90  132  116 

Y  s  47  46  53  31  64  43  75  40  80  55 


SOLUTION:  Hq:  p  =  0  ,  ^ :  p?  0 

(1)  Press  HD  ■  .05,  fr/sl,  jVs] . 

(2)  Enter  data  using  DEP  Sequence 

(3)  Press  [  A'  |  1  is  displayed  (reject  Hq) 

[x<t |  si  =  .0012735515 


9-  Two  types  of  electric  bulbs  are  observed  as  to  length  of  life,  with 


the  following  results:  n  =  46,  m  =  64,  x  =  1070  and  y  =  1041.  Are 


the  mean  lives  significantly  different?  (££  =  -10) 

§2HEI2»>  V  A Ay  ■  V  Ax*' Ay 

(1)  Press  |  s' |.  .1,  |a/S| .  0,  |b/S | . 

(2)  store  n  =  46  in  R^,  m  =  64  in  RQ^,  x  =  1070  in  R^  and 
y  =  1041  in  Ro8  . 

(3) ' [~E  |  0  is  displayed  (accept  Hq) 

|ScS.t  |  si  =  .5605970415 
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USER  GUIDE  FOR  FIVE  DISTRIBUTION  PROGRAMS 

INTRODUCTION:  The  purpose  of  these  five  programs  is  to  provide  accurate 
approximate  values  for  the  following  distributions:  Normal,  Binomial/ 
Multinomial,  Chi-square,  Student's  t  ,  and  F  . 

GENERAL  PROCEDURES:  Read  in  the  appropriate  program  and  proceed  as 
below . 

NORMAL  Distribution  -  Program  5 

Proceed  with  the  following  steps  in  any  order. 


Enter 

Press 

Display 

1. 

z 

A 

z)  =  normal  density  at  z 

2. 

z 

B 

P  (z)  =  P(Z  S  z) 

3* 

P(z) 

C 

Zp  =  normal  p^*"1  percentile 

4. 

z 

D 

Q(z)  =  P(Z  >  z) 

5. 

z 

E 

A(z)  =  P(ZS  |zj) 

BINOMIAL 

Distribution  -  program  6 

Perf 

om  the 

first  two 

steps  and  then  do  the  remaining  steps  in 

any 

order . 

Enter 

Press 

Display 

1  . 

n 

A 

n 

2. 

P 

3 

P 

3- 

k 

C 

f(k;n,p)  =  binomial  density  at  k 

4. 

k 

D 

F(k;n,p)  =  P(zSk) 

5. 

k 

E 

Q(k)  =  P(Z  >k) 

6. 

- 

A’ 

M 

7. 

- 

B’ 

a 
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MULTINOMIAL  Density  -  Program  6 


Perform  the 

following  steps  in  order. 

each  set  (n 

VVL)  for 

i—  1 1 2 ,  . . .  f  k .  (k  - 

Enter 

Press 

Display 

1.  N 

- 

2.  n^ 

£n^  (can 

pi 

[i7i] 

(Gan 

3. 

HD 

f (nl ,n2» *  * 

*  As  an  added  feature  the  following  capability  exists  in  Program  6. 
4.  Store  N  in  RQ^  and  press  |  S'  |  to  display  N.'  (N  —69). 

CHI -SQUARE  Distribution  -  Program  7 


Enter  Press  Display 

1-  X2  f  A  I  f(X2)  =  density  at  X2 

2.  x2  {jj  p(x2)  -  pcz^x2) 

3.  P(X2)  |  C  1  X2p(v)  =  ptn  percentile 

*  As  an  added  feature  the  Gamma  function  can  be  evaluated  as  below. 

4*  v  f  D~]  IW2) 

STUDENT'S  t  Distribution  -  Program  8 

Store  v  in  and  then  the  following  steps  in  any  order. 


Enter 

Press 

Display 

1. 

t 

CU 

f(t)  =  density  at  t 

2. 

t 

LlI 

P(t)  =  P(Z  s  t) 

3. 

P(t) 

m 

t  (v)  =  p^*1  percentile 

*  As  an  added  feature  the  Gamma  function  can  be  evaluated  as  below. 
5.  Store  v  in  and  press  [  D  |  to  display  J^(v/2) . 
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F  Distribution  -  Program  9 

Complete  Steps  1  and  2  prior 


Enter 

Press 

1. 

V1 

2. 

V2 

m 

3. 

> 

rH 

> 

LlJ 

4. 

P(F) 

CD 

PROGRAM  5 

LABELS  USED 

A 

BCD 

E 

PROGRAM  6 

LABELS  USED 

A 

BCD 

E 

PROGRAM  7 

LABELS  USED 

A 

BCD 

E 

PROGRAM  8 

LABELS  USED 

A 

3  G  D 

E 

PROGRAM  9 

LABELS  USED 

A 

BCD 

E 

PROGRAM  5 

REGISTERS  USED 

9  11 

25  26  29 

PROGRAM  6 

REGISTERS  USED  - 

■  ALL 

PROGRAM  7 

REGISTERS  USED 

01  09  10  11  13  14 

15 

PROGRAM  8 

REGISTERS  USED 

00  02  09  10  15  thru  24 
PROGRAM  9  REGISTERS  USED 


to  Steps  3  or  4  below. 
Display 

P(F)  =  P(Z  —  F) 

"tbi 

F  (Vi , v^)  =  p  percentile 

D’  E* 

A'  B*  G'  D' 

A'  B’  C'  D' 

B*  C’  E’ 

A'  B’  C  D* 

17  thru  23 


00  02  08  09  11  15  thru  29 


PROGRAM  5  ~  Sample  Problems 

1 .  Find  0(1 .960) . 

SOLUTION: 

Enter  I.96O  and  press {~~A~|  .  .0584409443  is  displayed 

2.  Find  P(1 .960)  =  P(Z  -  1 .960) 

SOLUTION: 

Enter  1 .960  and  press |  B  ],  .9750021748  is  displayed 

3.  Find  P(-l .960) . 

SOLUTION: 

Enter  -1 .960  and  press [  B  ),  .0249978252  is  displayed. 

4.  Find  z 

SOLUTION: 

Enter  .05  and  press j  G  |t  -1 .64521144  is  displayed. 

5.  Find  z  ^ 

SOLUTION: 

Enter  .95  and  pressj  G  j,  1.64521144  is  displayed. 

6.  Find  Q(l.282)  =  P(Z=*1.282) 

SOLUTION: 

Enter  1.282  and  pressj"~D  |,  .0999213886  is  displayed. 

7.  Find  A(l. 282)  =  P(ZS  (z|  ) 

SOLUTION: 

Enter  1.282  and  Pressj  i  |,  .8001572228  is  displayed. 


88 


H 


m 


PROGRAM  6  SAMPLE  PORBLEMS 

1.  Find.  f(k;n,p)  given  k  =  6,  n  =  10  and  p  =  .4  . 


SOLUTION! 

Enter  10 

press 

j  a  j  then  enter  .4  press  j"  B  j 

enter  6 

press 

|  G  |,  .111476?36  is  displayed. 

2. 

Find  P(8)  when  n 

=  10  and  p  =  .75  . 

SOLUTION: 

Enter  10 

press 

j  A  j  then  enter  .75  press  |  B  J 

enter  8 

press 

[  D  | ,  .7559747696  is  displayed. 

3- 

Find  Q(2)  when  n 

=  10  and  p  =  .25  • 

SOLUTION: 

Enter  10 

press 

|  A  j  then  enter  .25  press  J  B  j 

enter  2 

press 

|  E  j,  .474407196  is  displayed. 

4. 

Find  jj,  and  (J  for  a  Binomial  distribution  with  n  = 

SOLUTION: 

Enter  20 

press 

|  A  |  then  enter  .5  press  |~B  | 

press  |  A' 

]  M  = 

10  is  displayed,  then 

press  f B' 

](JZ 

=  2.236067977  is  displayed. 

5- 

Find  f  (n^) 

for  the  following  Multinomial  case. 

n±  1  2 

3 

4 

*3 

.2 

.1 

SOLUTION: 

Enter  10  press  [ID  then 

enter  n^  |x%t  j  p^  [r/s|  etc.  and 

press  HH]  ■  .00036288  is  displayed. 
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then 


then 


then 


20  and  p 


then 


PROGRAM  7  Sample  problems 


1.  Find  f(25)  where  v  =  15 . 

SOLUTION; 

Store  15  in  R^  ,  then  enter  25  and  press  m 
.0134298528  is  displayed. 

2.  Find  P(25)  where  v  =  15. 

SOLUTION i 

Store  15  in  R^  ,  then  enter  25  and  press  [  B  ) 
.9500565864  is  displayed. 

3.  FindX2<9g(l6). 

SOLUTION; 

Store  lo  in  R^  ,  then  enter  .99  and  press  [  G  | 
31.99987803  is  displayed. 

4 .  Find  IV)  • 

SOLUTION: 

Enter  16  press  |  D  |,  5040  is  displayed. 

PROGRAM  8  Sample  Problems 

1.  Find  f(2)  where  v  =  20 . 

SOLUTION; 

Store  20  in  R^  ,  then  enter  2  and  press  f  A  | 
.0580872152  is  displayed. 

2.  Find  P(-.860)  where  v  =  20 . 

SOLUTION; 

Store  20  in  R^  ,  then  enter  -.860  and  press  [  B  | 
.199990431  is  displayed. 


3.  Find  t<9Q(l4). 

SOLUTION! 

Store  14  in  R.,  ^  ,  then  enter  .90  and  press  m 

I. 345266653  is  displayed. 

4.  Find  P(4.5). 

SOLUTION! 

Store  9  in  1  then  press  m 

II. 6317284  is  displayed. 

PROGRAM  9  -  Sample  Problems 

1.  Find  P(2.52)  where  ^  =  5  and  v£  =  10. 

SOLUTION! 

Enter  5  and  press  [  x$t  [,  then  enter  10  and  press  |  A  ] 
enter  2.52  and  press  (  B  j,  .8998498088  is  displayed. 

2.  Find  P( -397)  where  =  10  and  v  2  =  5. 

SOLUTION! 

Enter  10  and  press  fx<t  j,  then  alter  5  and  press  |  A  | 
enter  .397  and  press  j  B  |,  .1002517964  is  displayed. 

3 .  Find  F. 975(1° ,30)  • 

SOLUTION! 

Enter  10  and  press  |  xVb  ),  then  enter  30  and  press  CD 
enter  .975  and  press  CD'  2.511201637  is  displayed. 


APPENDIX  B 


SELECTED  CHI-SQUARE  INVERSE  CDF  APPROXIMATIONS 


DEGREES 

OF 

FREEDOM 

TABLED 

VALUE* 

(—  TYPE 
APPROXIMATE 
PERCENTILE 

ACTUAL 

PROBABILITY 

1 —  TYPE  II  — 1 
APPROXIMATE  ACTUAL 

PERCENTILE  PROBABILITY 

X2  05  1 

.0039 

.0000 

.0001 

.0039 

.0499 

A.  .03  2 

.1026 

.0789 

.0387 

.1026 

.0500 

3 

.352 

.3280 

.0453 

■3531 

.0503 

4 

•  711 

.6900 

.0474 

.7115 

.0501 

5 

l.l  5 

1 .1277 

.0484 

1.1460 

.0500 

10 

3.94 

3.9307 

.0496 

3.9404 

.0500 

13 

7.26 

7.2543 

.0498 

7.2610 

.0500 

20 

10.83 

10.8455 

.0499 

10.8508 

.0500 

30 

18 .49 

18.4885 

.0499 

18 .4927 

.0500 

40 

26.51 

26.5056 

.0499 

same  as 

TYPE  I 

60 

43.19 

43.1843 

.0500 

It  11 

ll 

120 

95-70 

95.6999 

.0500 

i«  ll 

It 

X2  90  1 

2.71 

2.6395 

.8958 

2.7067 

.9001 

4.61 

4.5596 

.8977 

4.6052 

.9000 

3 

6.25 

6.2146 

.8984 

6.2517 

.9000 

4 

7-78 

7.7480 

.8987 

7.7797 

.9000 

5 

9.24 

9.2086 

.8990 

9.2366 

.9000 

10 

15-99 

15.9688 

.8995 

15-9873 

.9000 

15 

22.31 

22.2930 

•  8997 

22.3072 

.9000 

20 

28.41 

28.4003 

.8997 

28 .4120 

.9000 

30 

40.26 

40.2473 

.8998 

40.2560 

.9000 

40 

51.81 

51 -7981 

.8999 

same  as 

TYPE  I 

60 

74.40 

74.3923 

.8 999 

If  ll 

II 

120 

140.23 

140.2309 

.9000 

ll  ll 

ll 

Y 2  on  c  1 

7.88 

7.9071 

•  9951 

7.8815 

•  9950 

A  >995  2 

10.60 

10.6753 

•9952 

10.5966 

■9950 

3 

12.84 

12.9227 

•9952 

12.8377 

■  9950 

4 

14.86 

14.9437 

•  9952 

14.8602 

•9950 

5 

16.75 

16.8303 

.9952 

16.7497 

.9950 

10 

25.19 

25.2557 

•  9951 

25-1884 

.9950 

15 

32.80 

32.8604 

•9951 

32.8015 

.9950 

20 

40.00 

40.0502 

•9951 

39.9970 

•9950 

30 

53.67 

53- 71 81 

•  9951 

53.6721 

•9950 

40 

66.77 

66.8076 

•  9950 

same  as 

TYPE  I 

60 

91.95 

91 .9879 

•  9950 

it  ll 

It 

120 

163-64 

163.6777 

•  9950 

11  11 

11 

*  [Ref.  2:  p.  46 5] 
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SELECTED  STUDENT'S  t  INVERSE  CDF  APPROXIMATIONS 
DEGRESS  I —  TYPE  I - 1  j—  TYPE  II  — ) 


OF  TABLED  APPROXIMATE  ACTUAL  APPROXIMATE  ACTUAL 

FREEDOM  VALUE*  PERCENTILE  PROBABILITY  PERCENTILE  PROBABILITY 


*.60  1 

.325 

.2568 

.5800 

.3249 

.6000 

2 

.289 

•  2593 

•  5902 

.2869 

•  5994 

3 

.277 

.2556 

.5926 

•  2754 

.5996 

4 

.271 

.2547 

•  5942 

.2699 

•5997 

5 

.267 

.2542 

•  5953 

.2667 

•  5999 

10 

.260 

•2535 

•  5975 

.2660 

•  5999 

15 

.258 

.2533 

•  5983 

.2578 

.6000 

20 

<2  57 

.2532 

•  5987 

.2567 

.6000 

30 

.256 

.2531 

•  5990 

.2556 

.6000 

40 

•255 

.2531 

•  5992 

.2550 

.6000 

60 

.254 

•  2530 

•  5994 

.2545 

.6000 

120 

.254 

.2530 

•  5996 

.2539 

.6000 

^  .90  * 

3.078 

1.7878 

.8377 

3-0777 

.9000 

y  2 

1.886 

2.1065 

•9151 

1.8764 

•8993 

3 

I.638 

1 .6292 

.8991 

I.6383 

.9001 

4 

1.533 

1.5092 

.8971 

1.5347 

.9002 

5 

1.476 

1 .4420 

.8969 

1 .4772 

.9002 

10 

1.372 

1.3581 

.8979 

1.3726 

.9001 

15 

1.341 

1-3311 

•  8985 

1.3408 

.9000 

20 

1.325 

1.3182 

.8988 

1.3255 

.9000 

30 

1.310 

1.3057 

.8992 

1.3105 

.9000 

40 

1.303 

1 .2996 

.8994 

1.3031 

.9000 

60 

I.296 

1 .2936 

.8996 

1.2958 

.9000 

120 

1.289 

1 .2876 

.8998 

1.2886 

.9000 

^  .995  ^ 

63.657 

6.6860 

•  9527 

63.6567 

•  9950 

9.925 

9.2732 

•  9943 

9.7472 

•  9948 

3 

5-841 

5- 3978 

•  9938 

5.8223 

.9950 

4 

4.604 

4.4230 

.9943 

4.5990 

•  9950 

5 

4.032 

3.9585 

.9946 

4.0299 

.9950 

10 

3.169 

3.1961 

•  9952 

3.1702 

•  9950 

15 

2.947 

2.9771 

•  9953 

2.9476 

•9950 

20 

2.845 

2.8726 

•  9953 

2.8460 

•  9950 

30 

2.750 

2.7712 

•  9953 

2.7504 

•  9950 

40 

2.704 

2.7215 

•  9952 

2.7047 

•  9950 

60 

2.660 

2.6724 

■  9952 

2.6604 

•  9950 

120 

2.617 

2.6240 

•  9951 

2.6175 

•9950 

*  [Ref.  2:  p.  464] 
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SKT.-RGTSP  F  INVERSE  CDF  APPROXIMATIONS 


DEGREES 

OF 

freedom 

TABLED 

VALUE* 

/ - TYPE 

APPROXIMATE 

PERCENTILE 

I - » 

ACTUAL 

PROBABILITY 

l —  TYPE 
APPROXIMATE 
PERCENTILE 

II  —”i 

ACTUAL 

PROBABILITY 

vl  v2 

F  1  i 

.0062 

.00619 

.0500 

.05  2 

.0050 

- 

- 

.00392 

.0443 

5 

.0043 

- 

- 

.00391 

.0475 

15 

.0041 

- 

- 

.00391 

.0490 

30 

.0040 

- 

- 

.00391 

.0495 

120 

.0039 

.00391 

.0498 

F  rsc  5  1 

•  151 

_ 

_ 

.1515 

.0501 

.05  2 

.173 

.1340 

•0315 

.1611 

.0442 

5 

.198 

.1949 

.0495 

.1980 

.0500 

15 

.216 

.2110 

.04?4 

.2165 

.0500 

30 

.222 

•  2157 

.0469 

.2223 

.0500 

120 

.227 

•  2195 

.0464 

.2272 

.0499 

F  n, i5  1 

.220 

_ 

_ 

.2190 

.0495 

•°5  2 

.272 

.2054 

.0235 

.2435 

.0378 

5 

•  34  3 

.3432 

.0489 

•3446 

.0500 

13 

.416 

.41 57 

.0498 

.4161 

.0500 

30 

.445 

.4442 

.0496 

.4451 

.0500 

120 

.473 

.4713 

.0492 

.4730 

.0500 

F  n ,  30  1 

.240 

- 

.2388 

.0496 

.05  2 

.302 

.2260 

.0207 

.2671 

.035^ 

5 

•  395 

•  3909 

.0483 

.3945 

.0499 

15 

.496 

.4962 

.0499 

.4963 

.0500 

30 

.543 

.5431 

.0499 

.5432 

.0500 

120 

•  594 

•  593  5 

.0497 

.59^0 

.0500 

f.o5  120  l 

•255 

- 

- 

.2547 

.0498 

.326 

.2421 

.0184 

.2855 

.0332 

5 

.437 

.4310 

.0474 

.4364 

.0498 

15 

•  571 

•  5705 

.0496 

•  5713 

.0500 

30 

.644 

.6431 

.0499 

.6434 

.0500 

120 

.740 

•  7397 

.0500 

•  7397 

.0500 

*  [Ref.  2:  pp.  472-485] 


SELECTED  F  INVERSE  CDF  APPROXIMATIONS 


DEGREES 

OF 

FREEDOM 

TABLED 

VALUE*' 

1 —  TYPE 
APPROXIMATE 
PERCENTILE 

I  — l 

ACTUAL 

PROBABILITY 

1 —  TYPE 
APPROXIMATE 
PERCENTILE 

II— 1 

ACTUAL 

PROBABILITY 

V1  V2 
F.995  1  1 

16200 

16210.7 

.9950 

2 

198 

- 

- 

135.65 

.9925 

5 

22.8 

- 

- 

21 .125 

.9941 

15 

10.8 

- 

- 

11 .063 

•  9954 

30 

9.18 

- 

- 

9.361 

.9954 

120 

8.18 

“ 

8.233 

•  9951 

F-995  5  1 

23100 

- 

- 

25615.1 

•  9953 

2 

199 

554.13 

.9982 

276.56 

•  9964 

C, 

> 

14 .9 

15.872 

.9956 

15.008 

.9951 

15 

5.3? 

5.409 

•9951 

5-373 

•9950 

30 

4.23 

4.245 

•9951 

4.228 

•9950 

120 

3.55 

3-558 

■9951 

3.549 

•  9950 

7  1  S  1 

*  .995  °  g 

24600 

- 

- 

25615.2 

•9951 

199 

572.87 

.9983 

285.06 

•  9965 

5 

13-1 

14.009 

■9957 

13.228 

•  9951 

15 

4.07 

4.082 

•  9951 

4.070 

.9950 

30 

3.01 

3-007 

.9950 

3.006 

■  9950 

120 

2.37 

2.372 

•9950 

2.373 

■9950 

F.995  30  l 

25000 

- 

- 

25615.2 

.9951 

199 

578.74 

•9983 

287-53 

.9965 

e, 

y 

12.7 

13.530 

■9951 

12.747 

.9951 

15 

3-69 

3.700 

.9951 

3-687 

•9950 

30 

2.63 

2.629 

•  9950 

2.628 

•9950 

120 

1.98 

1 .984 

•  9950 

1 .984 

■9950 

F.995120  3 

25400 

_ 

_ 

25615.2 

■  9950 

199 

583.29 

•  9983 

289 .43 

.9966 

3 

12.3 

13.159 

•9958 

12.373 

•  9951 

15 

3-37 

3-387 

•  9951 

3-373 

■  9950 

30 

2.30 

2.302 

•9950 

2.300 

■  9950 

120 

1.61 

1 .606 

•  9950 

1 .606 

■9950 

COMPUTER  LISTINGS 


PROGRAM  1  ONE-POPULATION  CONFIDENCE  INTERVALS 
LABEL  ADDRESSES  non  d'O  c t n 


001  71  SBR 

046  14  D 

0  59  12  B 

114  16  fl ' 

1  O  1  p  p 
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1  P,n  33  n- 


o  *  n  J.  ! 
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450  1 3  C ' 

464  61  GTD 

480  43  RCL 


PROGRAM  LISTING 


03  03 

43  RCL 
0  1  0  1 


020 

42  STD 

063 

li- 

I N  V 

021 

21  21 

064 

1-, 

EQ 

022 

7 1  S  B  R 

0  6  5 

00 

0  0 

023 

A  •"«  ri  « 

“*t  ■  •=  L-  L, 

0  6  6 

69 

0  9 

024 

42  STD 

067 

71 

S  B  R 

025 

04  04 

! J  0  S 

71 

SBR 

026 

43  RCL 

1 J  0 3 

53 

027 

“i  4  -i  4 

ii  ^1 

070 

.J  .‘l 

f’ 

028 

!-( p;  4- 

071 

H  p  ~- 

.4  ~! 
b+  “« 

01 

o  .j 

RCL 

M  1 

0  y 

U  3  0 

54  > 

0  7  3 

'  J  i 

U  3 1 

■1  :~i  !~  «. »“• 

•-+  L- 

074 

01 

032 

2  0  2  0 

075 

54 

0  3  3 

75 

0  7  6 

65 

034 

Ij  £  5^ 

077 

43 

F“r  I-*  1 

K  i_.  L 

035 

54  > 

078 

i  i 

1  1 

0  3 

42  STD 

079 

n  j 

- 

0  3  7 

•~i  -i  •  J  -j 

080 

c* 

f 

71  SBR 

081 

24 

“•  r- 

L.-  L. 

0  3  9 

4  3  R  C  L 

n  8  p 

+ 

040 

43  RCL 

0  3  3 

53 

041 

0  4  0  4 

084 

.4  -i 

HO 

O  *•  * 
r.  «_•  l. 

042 

42  STO 

085 

03 

0  3 

i  j  4 

■»  4  4  4 

i  1  i  i 

0  *3  6 

75 

- 

044 

92  RTN 

U  8  7 

h*  •— * 

RCL 

045 

76  LBL 

0  8  8 

01 

0 1 

046 

14  n 

0  8  9 

54 

j 

047 

25  CLR 

090 

5  4 

) 

048 

04  4 

091 

54 

0  4  9 

f,  q  n  p 

092 

*”i 

>  T 

050 

H  —?  4 

if  if 

0  9  3 

C.  -~t 

< 

U  5  1 

2  9  C  P 

0  9  4 

.4  ~i 
H*  .  j 

r>  “• ; 

052 

47  CMS 

0  8  0 

fi  :-i 

053 

91  R/S 

0  9  6 

l'  5 

- 

054 

78  2  + 

097 

*■+  o 

R  C  L 

055 

61  GTD 

098 

01 

01 

056 

00  00 

099 

85 

+ 

0  5  7 

53  53 

100 

0  1 

i 

i 

058 

76  LBL 

101 

54 

059 

12  B 

1 02 

65 

060 

29  CP 

103 

.•!  •*i 

“f  mZ* 

R  U  L 

0  6 1 

43  RCL 

104 

1  3 

i  J 

U  £  3 

11  11 

105 

55 

- 

PROGRAM  2  TOO -POPULATION  CONFIDENCE  INTERVALS 

026  61  GTD 

LABEL  ADDRESSES  027  79  ~ 

028  42  STD 

UUi  19  I 029  03  08 

U  £  b  61  G  T  0  n  n  o  t  kj 


PROGRAM  LISTING 

0  0  0  7  6  L  B  L 

001  19  D ' 

002  43  RCL 

003  01  01 

004  42  STD 

005  13  13 

006  43  RCL 

007  02  02 

003  42  STD 

8  09  14  14 

010  43  RCL 

011  03  03 

012  42  STD 


42  STD 
1 3  1 3 

22  INV 
79  x 

42  STD 
17  17 

61  GTD 
0  4  0  4 

6  2  6  2 
76  LBL 


026 

61  GTD 

074 

42 

STD 

027 

79  X 

075 

19 

19 

023 

42  STD 

076 

54 

*! 

029 

08  08 

077 

34 

r  *..* 

4 

030 

£2  INV 

073 

42 

STD 

031 

79  x 

079 

•”i  o 

2 

032 

“i  ■”«  .* 

."i  .Ji  >- 

030 

43 

RCL 

033 

42  STD 

031 

19 

19 

034 

07  07 

032 

92 

RTH 

035 

43  RCL 

083 

7  6 

LBL 

036 

1  3  1 3 

084 

95 

= 

037 

42  STD 

085 

A  ~i 

HO 

RCL 

033 

04  04 

086 

i  ■? 

12 

039 

43  RCL 

037 

hi  5 

040 

14  14 

088 

43 

r~\  I 

041 

42  STD 

0  8  9 

1 1 

4  ^ 

1  i 

042 

05  05 

090 

O  c; 

+ 

043 

43  RCL 

091 

43 

RCL 

044 

15  15 

092 

10 

10 

045 

75 

093 

5  4 

’j 

0  4  6 

01  1 

094 

92 

RTH 

047 

54  > 

095 

7*  S 

LBL 

043 

65  x 

096 

13 

j“* 

049 

43  RCL 

097 

£9 

CP 

050 

I  T  IT 

II  i  ! 

093 

43 

RCL 

051 

35  + 

099 

10 

10 

052 

53  < 

1 66 

67 

EQ 

053 

53  < 

1 0 1 

01 

01 

054 

43  RCL 

102 

•n*  ■? 

055 

03  03 

1  03 

43 

RCL 

0  5  6 

n- 

(  .j 

1  04 

15 

15 

057 

01  1 

105 

67 

EQ 

0  5  3 

54  > 

106 

01 

01 

059 

6  5  x 

107 

50 

50 

u  t*  u 

43  RCL 

103 

43 

R  C  L 

061 

0  7  0  7 

109 

16 

i  i 

i  V 

062 

er  a 
•  J  H 

1 1 0 

67 

EQ 

063 

54  ) 

1  1 1 

01 

01 

064 

55  - 

112 

35 

35 

065 

53  < 

4  4  C' 

1  i 

85 

+ 

0  6  6 

43  RCL 

114 

43 

RCL 

0  67 

1 5  1  5 

115 

26 

063 

35  + 

116 

54 

i 

069 

43  RCL 

117 

5  5 

- 

070 

03  03 

1  t  0 

43 

RCL 

07 1 

75  - 

119 

19 

19 

!J  )■'  ^ 

02  £ 

1  20 

54 

) 

073 

54  > 

121 

34 

n< 

PROGRAM  2  Continued 


26b 

29 

29 

314 

n  pt 

.“1 

362 

02 

M  d 

267 

85 

+ 

315 

65 

3  6  3 

55 

- 

2  6  S 

93 

316 

43 

RCL 

364 

43 

RCL 

969 

no 

n 

317 

29 

29 

365 

04 

04 

270 

01 

l 

313 

45 

Y 

3  6  6 

54 

j 

271 

00 

0 

319 

03 

3 

367 

42 

STD 

“i  -?  •-* 

03 

320 

54 

j 

3  6  8 

0  8 

08 

•-»  "7  -•j 

02 

321 

54 

y 

369 

94 

+  /  - 

*2.  l*'  *4 

03 

Cj 

-»  0  0 

94 

+ /  “ 

370 

85 

+ 

Qi  *7  C| 

65 

3  2  3 

85 

+ 

3  7 1 

01 

1 

276 

43 

RCL 

324 

43 

RCL 

3  «■'  2 

54 

) 

9 

yQ 

325 

29 

99 

•2*  i  -2* 

65 

X 

•n  -7  i^i 

2 

X  ^ 

326 

54 

j 

374 

A  -“I 
*-r  .  “1 

RCL 

279 

5  4 

•-»  't* 

i1  d  r 

42 

STD 

3  7  5 

Q9 

i  IS 

2  ?:n 

55 

4. 

1 1 

1 1 

3  7  6 

55 

- 

O  O  J 

53 

( 

3  2  9 

92 

RTN 

•2«  { 

43 

RCL 

i-»  “i 

01 

330 

7  6 

LBL 

-2'  {  O 

04 

04 

•n  O  "■ 

35 

-1  •-«  i 

1  '!’ 

1  i— 

hJ 

379 

54 

y 

294 

n  i 

>} 

I1  *2>  O 

99 

3  8  0 

35 

+ 

•-«  -i  ir 

21  O  _l 

93 

43 

RCL 

381 

43 

RCL 

.-i  .  ~ 

.  ’  ”i  i  i 

fid 

4 

334 

1  1 

4  4 

i  i 

■j  O  •”! 

05 

05 

no 

335 

9  9 

IHV 

l“l  'Ij 

54 

y 

w  '  !_* 

U  d 

336 

67 

EQ 

384 

•"1  4 

•2*  "t 

rx 

2  9  ^ 

07 

1 

2’  '2* 

03 

03 

385 

42 

STD 

290 

0  3 

O 

■?  O  O 

40 

40 

386 

1  *2* 

X  u- 

12 

9  9  1 

C; 

2j_3_9 

17 

B  * 

•-t  ;“»  ”7 

•_=  i 

43 

RCL 

292 

65 

340 

43 

RCL 

■j  O  O 

18 

18 

2  9  3 

A  •*> 

‘t  2» 

RCL 

341 

01 

01 

389 

75 

~ 

294 

99 

29 

3  42 

55 

-i- 

3  9  0 

485 

RCL 

295 

35 

2*  2' 

43 

RCL 

391 

08 

08 

296 

•I:  ■? 

. 

”t  9  A 

- f  *-? 

03 

03 

392 

54 

> 

297 

0 1 

i 

345 

54 

y 

393 

42 

STD 

2  9  8 

0  3 

346 

42 

STD 

394 

10 

1  0 

2  Q  Q 

135 

3 

“»  .«{  "7 

2'  4  | 

4 

i  O 

4  <~t 

1  0 

3  9  5 

4  £ 

i 

H  * 

2  0  0 

f  1 9 

-1 

**;  4 !"! 

94 

+ - 

0  9  6 

92 

RTN 

3  n  i 

06 

6 

349 

0  0 

+ 

3  9  7 

7  6 

LBL 

302 

119 

9 

350 

01 

4 

i 

3  9  8 

10 

c  ■ 

303 

65 

2;  R  1 

54 

> 

3  9  9 

43 

RCL 

304 

43 

RCL 

352 

65 

X 

400 

08 

0  8 

305 

2_3 

29 

“t  C-  •-( 

2*  J  -2' 

43 

RCL 

401 

42 

STD 

2.  n  9 

■? 

v  2. 

854 

18 

1  8 

402 

17 

307 

85 

4- 

355 

55 

-1- 

403 

43 

RCL 

308 

93 

356 

43 

RCL 

404 

i  O 

13 

309 

0  n 

n 

*■»  r  -7 

2»  -J  r 

03 

03 

405 

42 

STD 

310 
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0 

358 

54 

’) 
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07 

■  J  i 

311 

01 

1 
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42 

STD 

407 

43 
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05 
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no 

0 
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43 
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29 
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410 

43 

RCL 
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03  3 
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29 

CP 
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459 
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